Exact prefactors in static and dynamic correlation 
functions of ID quantum integrable models: 
applications to the Calogero-Sutherland, 
Lieb-Liniger and XXZ models 

17 October 2011 

Aditya Shashi 

Department of Physics and Astronomy, Rice University, Houston, Texas 77005, 
USA 

Milosz Panfil 

Institute for Theoretical Physics, Universiteit van Amsterdam, 1090 GL 
Amsterdam, The Netherlands 

Jean-Sebastien Caux 

Institute for Theoretical Physics, Universiteit van Amsterdam, 1090 GL 
Amsterdam, The Netherlands 

Adilet Imambekov 

Department of Physics and Astronomy, Rice University, Houston, Texas 77005, 
USA 

Abstract. In this article we demonstrate a recently developed technique which 
addresses the problem of obtaining non-universal prefactors of the correlation 
functions of ID systems at zero temperature. Our approach combines the effective 
field theory description of generic ID quantum liquids with the finite size scaling 
of form factors (matrix elements) which are obtained using microscopic techniques 
developed in the context of integrable models. We thus establish exact analytic 
forms for the prefactors of the long-distance behavior of equal time correlation 
functions as well as prefactors of singularities of dynamic response functions. 
In this article our focus is on three specific integrable models: the Calogero- 
Sutherland, Lieb-Liniger, and XXZ models. 
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1. Introduction 

One finds for one-dimensional (ID) quantum mechanical systems of interacting bosons, 
fermions and spins, that the highly restrictive kinematics leads to many examples of 
exactly solvable models. Consequently there often exist explicit expressions for the 
wavefunctions of such models, and for some of their bulk properties like energy and 
the excitation spectrum. However in the context of condensed matter physics one 
is interested in the correlation functions of these systems, which characterize their 
response to external probes. Unfortunately the complexity of the wavefunction makes 
the problem of calculating correlation functions notoriously difficult. 

On the other hand, the effective low energy physics of a broad set of ID 
quantum liquids can be described using a hydrodynamic approach known as the 
Luttinger liquid theory [TJ [2j [3j 2J [5] . This theory predicts the universal long-range 
behavior of the correlation functions of such systems. One obtains, for example, 
asymptotic expressions for equal-time correlation functions at zero temperature as 
series expansions with power laws controlled by a dimensionless Luttinger liquid 
parameter K > 0, see Eqs. While the "universal" parameter K is 

related to thermodynamic properties and can be easily extracted from numerical or 
exact solutions, the "non-universal" (model-specific) prefactors in the series expansion, 
e.g. A m , B m ,C m , D m , see Eqs. ([I])-©, ©-©, are usually not known analytically 
except for a few cases [HI 13 [HI 13 HH1 HH H31 H21 EH- These prefactors set the actual 
scale of observable correlations, and are an impediment to establishing the complete 
asymptotic behavior of the static correlation functions of a given ID quantum model. 

In a recent publication |15) . by combining the analysis of the Luttinger liquid 
Hamiltonian with the finite-size properties of certain matrix elements (form factors), 
a general technique for calculating these non-universal prefactors was developed for 
a generic ID quantum liquid. Moreover, it has been shown recently 16j-|34] that 
dynamic response functions generically have singularities which can be described by 
effective Hamiltonians of impurities moving in Luttinger liquids. Analysis of the finite- 
size properties of these effective Hamiltonians can also be used to obtain prefactors of 
various dynamic response functions. 

In what follows, we apply this technique to the calculation of prefactors of the 
correlation functions of integrable models, focusing specifically on three integrable 
models: the Calogero-Sutherland model (CSM) [35] [36] [37] of fermions interacting 
via a long range inverse-squared-distance potential, the Lieb-Liniger model [38] of 
bosons with pairwise contact interactions and the XXZ model [39 , 40J of a Heisenberg 
spin chain. These models have the advantage of being solvable [37J [351 Ell HI] and 
additionally have readily available expressions for finite size form factors gSJ HSJ 
[47l [48] . which we will investigate to obtain analytic expressions for prefactors of 
their correlation functions valid in the thermodynamic limit. Furthermore, the Lieb- 
Liniger model has been realized with ultracold atomic gases |50j . and its correlation 
functions can be measured using interference (STJ [52j [53] , analysis of particle losses [54] . 
photoassociation [S3], or Bragg and photoemission spectroscopy 55], and density 
fluctuation statistics [56 . The XXZ model has enjoyed a long history as a well 
studied model in statistical physics (43] [57] [58] . 

This article is organized as follows. Sec. 2 we outline the analysis of the effective 
field theory description of ID quantum systems. In Sec. 2.1, we use the linear 
Luttinger liquid theory to work out the connection between prefactors of equal-time 
correlation functions and lowest energy form factors. In Sec. 2.2, we show that 
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analysis of the effective field theory of impurities moving in Luttinger liquids gives 
the relations between form factors and prefactors of dynamic response functions. In 
Sec. 3 we present details of the calculation of the prefactor in the density structure 
factor (DSF) in the vicinity of the threshold singularity for the CSM, by working out 
the thermodynamic limit of form factors. This is to serve as an introduction to the 
more technically involved calculations presented in Sec. 4 where we obtain several 
prefactors of static and dynamic correlation functions of the Lieb-Liniger model, with 
details on the calculation of thermodynamic limits of form factors of this model. In 
Sec. 5 we present the derivation of the prefactors of the equal time spin-spin correlator 
of the XXZ spin chain with a finite magnetization. Additional technical details are 
contained in the Appendixes. 



2. Results from effective field theory 

The Luttinger liquid theory [TJ [3J 01 [SJ |3] predicts the behavior of the correlation 
functions for spinless bosons and fermions of density p$, when pox > 1 as (here 
k F = 7rpo) 



(p(x)p(O)) K An cos(2mk F x) 

Pi * 2(w) 2+ 4' 1 (p xf m2K ' U 

(^ B (x)tp B {0)) ^ \ - B m cos(2mk F x) 

n n ~ 2-^ I N2m 2 K+l/(2if) ' W 

Po m >o(PQ x ) 
(j> F (x)$ F (Q)) ^ C m sin [(2m + l)k F x] 

Here p is the density operator, V'b(^'f) is the bosonic (fermionic) annihilation operator. 
The Hamiltonian describing these correlations is written as (we follow notations of 
Ref. [5]) 

Ho = iJ dx { K(V6)2 + i, (V0)2 ) ' (4) 

where v is the sound velocity, the canonically conjugate fields (f>(x),6(x) have the 
commutation relation [4>(x), V0(a;')] = inS(x—x / ), and the components of the fermionic 
(bosonic) fields with momenta (2m + 1/2 ± l/2)k F are written as 

lpF(B)(x,t) ~ e i(2m+l/2±l/2)[fe F x-0(x,i)]+ie(x,t) ) (g) 

where each component is defined up to a non-universal prefactor. For repulsive bosons, 
one has K > 1, while for repulsive (attractive) fermions K < 1(> 1). In Eqs. ([T])- 
(|3"T). we included only slowest decaying power laws for each oscillating component. In 
principle, irrelevant corrections ( see e.g. Ref. [61]) to Hamiltonian (j4|) and operators 
([5]) generate various faster decaying power-law terms for each oscillating component 
in addition to the one presented above |62j . 

In the case of the spin chain, we consider specifically the XXZ model given by 
the Hamiltonian 

N 

H XXZ = Y, J (S*S* +1 + S?S? +1 + A5?5J +1 - hS*), (6) 

3=1 
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where Sj x = c"/2 is a spin half operator acting at site j, satisfying the algebra 
[S a ,S p ] = ie a p 7 S~<, and the spins at different sites commute. A is a measure of 
the anisotropy of the coupling in the z— direction and h is an external magnetic field 
aligned in the z— direction. One typically performs a Jordan- Wigner transformation 
to map the Hamiltonian ([6]) to that of a ID lattice of fermions with nearest neighbor- 
interactions. We can then calculate the asymptotic behavior of the correlators of the 
XXZ model, by first representing the spins in terms of fermions and then carry out the 
usual bosonization procedure |5j 3, 59] for the fermionic fields. A finite magnetization, 
s z = (S z ) is easily accommodated in this framework and enters as a chemical potential 
for the fermions, see e.g. Ref. [5]. Following this prescription one obtains the following 
bosonic representations for the spin operators: 

S Z (x,t) ~ S z - — + e l2 m[(s, + l/2)7r a; -0( £C ,t)]^ 
7T 

S + {x t) ~ e -i2m[(a,+l/2)x-^(x,i)]-»fl(a!,i) (y\ 

where again each oscillating component is defined up to a non-universal prefactor. 

Similarly to Eqs. (HJ-©, Hamiltonian Q then results in the following long 
distance behavior 



m>l 



IQ+I \q-(c\\\ E m cos(2m{s z + 1/2)ttx) 

(S + (X)S (0)) = ^ X 2m-K+1/(2K) ■ ( 9 ) 

m>0 

We note however, that in the presence of the lattice, the definition of the prefactors 
for higher m requires more care. Indeed, for rational fillings 1/2 + s z = p/q, the 
momenta of oscillating terms for m% = m2(mod q) are not distinguishable, since 
momentum is no longer a good quantum number, but quasimomentum is. Thus in 
principle "subleading" terms from smaller mi need to be distinguished from the leading 
term for mi. While this might not be possible for rational fillings (in particular, half- 
filling), and some rational values of the Luttinger liquid parameters K, it should be 
possible for generic filling and interaction strength. In what follows, we will focus 
on the case of generic interaction strength and finite magnetic field to avoid the 
complications related to half-filling (the case of half-filling for the XXZ model was 
considered in Ref. [T3] by different techniques). 

2.1. Prefactors of equal-time correlators from the Luttinger liquid theory 

Let us consider a system of interacting bosons. Using the resolution of the identity in 
the expectation value (tp^ B (x,t)ipB(0)), we get 

k.uj 

where (k,w\ipB\N) is a form factor of the annihilation operator, \k,u>) denotes an 
eigenstate of N — 1 particles with momentum k and energy u>, and \N) is the ground 
state of N particles. For a finite system, k and w are not continuous, but will be 
quantized and consequently the spectral function is a collection of delta functions 
in (k,u>). We will now obtain a similar representation from the Luttinger liquid 
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theory and compare it with Eq. (|10l) to obtain the non- universal prefactors B m . 
Hamiltonian (jj) can be written using left- and right-moving components (puR) = 
OvK^tp/^/K, [4] which dictates the time dependence of the cos (2mkFx) component 
of {^B{x,t)tpB{^, Q))/po, at p \x ±vt\ > 1 as 

' cos (2mkFX) 



-l) m B m p ' 



(11) 



(i(vt + x)+0f L (i{vt-x)+0f R ' 
where Hum = m 2 K±m+l/4K > 0. The coefficients B m appeared in Eq. (fTTj) because 
we relate the t — limit of (ip B (x,t)tpB{0,0)) / po to the right hand side of Eq. ([2]). 
The two factors in the denominator describe contributions from left (right)-going 
excitations which propagate with velocities =pu, and signs of the infinitesimal shifts 
in the denominators ensure that only excitations with negative (positive) momenta 
can be created at the respective branches. For a finite system with periodic boundary 
conditions on a circle of length L, conformal invariance dictates (see e.g. Ref. [4]) , 
the right-going component (i(vt — x) + 0) fJ ' R is replaced by 



a i7r(vt— x) I L 



iL sin 



7r(vt — x) 



^2 C{n r ,p R ) 



n r >0 



^2i7rn r (x — vt) / L 

(i/27r)^H ' 



C(n r ,p R ) 



(12) 



(13) 



T(p B )T(n r + 1)' 

and similarly for left- going components with p Rl n r substituted by pL,ni- These 
equations lead to nontrivial predictions for the exact scaling of the form factors of 
a model describing interacting ID bosons, e.g. the Lieb-Liniger model. By comparing 
the finite size excitation spectrum of Hq with the exact solution of the Lieb-Liniger 
model, up to oc 1/L terms, we identify = F27rn/( r )/L as the total momenta of excitations 
created near the left (right) quasi-Fermi points, i.e. the edges of the distribution of 
quasimomenta characterizing an eigenstate of the model. Considering n r = ni = 
then leads to the scaling law 



;m,JV-l|^B|0,JV> 



= {-l) m B mPo ( 2ir V 
2-<5o,m \PqLJ 



(14) 



where \m,N) denotes an eigenstate of N bosons having center of mass momentum 
2mkp,. We see that as a consequence of the criticality of the Luttinger liquid, form 
factors of the annihilation operator have nontrivial scaling with the system size, and 
the prefactors of these nontrivial powers of L are directly related to the prefactors of 
the correlation functions. 

For density correlations, field correlation functions for fermions and correlators 
of spins, similar relations can be worked out as long as we represent the relevant 
operators in the bosonized language, and are given by 



[m,N- 1\4>f\N) 



\(m,N\p\N)\' 



C m po ( 2ir 
2(-l) m \p L 



2K 



(m,N\S z \N) 



A m p ( 2tt \ 
2 \p Lj 

2m 2 K 



2m 2 K 



El 

2 



(15) 



(16) 



(17) 
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2 f_iynR /9 _\ 2m 2 AT+l/(2X) 

(m,iV -l|5-|iV)| » I ^ ^ j . (18) 

Eqs. (fI3 |) -(fj"8" ]l allow one to evaluate the prefactors in Eqs. ©-©, by 

identifying a single, simplest "parent" form factor for each of the operators p, ipB, tpF, 
S^andS-. 

We note here that it recently came to our attention that a method to obtain 
prefactors of equal-time correlators, similar in spirit to ours is presented in Ref. [64], 
however recent theoretical developments allow us to push the technique further and 
obtain e.g. the precise splitting of spectral weight among low energy form factors, and 
also to obtain prefactors of dynamic response functions. 

Field theoretical considerations allow to fix not only the form factor with n r = 
ni = 0, but also the form factors for all low-energy states. However, since for n r (^ > 1 
states can be degenerate, one needs to understand how the spectral weight C(n r , pn) 
is split between different form factors. This question can be answered by matching 
contributions from each form factor with the free fermionic quasiparticle representation 
of the Luttinger Liquid [35]. Such representation has the same degeneracies as the 
exact solution, and we calculated [TS] its form factors using the results of Ref. [63 . 
For a state with k particle-hole excitations near the right quasi-Fermi point specified 
by integers p\ > ... > Pk > (particles) and qi < ... < qk < (holes), with total 
momentum (2tt/L) Y]. (pj — q^ — 2i:n r / 'L, we obtain that the ratio of its form factor 
to the one with n r — ni = equals 

/({Pi.ft}) = Dci, r , (—t—)l\f + (pi)r(qi), (19) 
\Pi Id. 



i<k 



where 



f iP) = Tr, ; — vFZ 7T\ > f W = 



r(-v/^)r(p + i) !J v ^ r(i + VM5)r(-g) 

Normalization of the spectral weight leads to the following "multiplet summation rule" 
(see Ref. [15] and appendix of the same for details): 

E |/(fe,%})| 2 = C(n r ,/ifl). (20) 

When n\ ^ 0, contributions from the left quasi-Fermi point are accounted for similarly, 
and the total form factor is a product of these two terms. 



2.2. Prefactors of singularities in DRFs using three subband model 

We now apply the techniques described above to the prefactors of singularities in 
dynamic response functions [15]. For more comprehensive discussions of singularities 
in the response functions of ID quantum liquids as well as the field theoretical 
description which captures these phenomena, see Refs. [16]- [34]. For concreteness we 
focus here on the Lieb-Liniger model of bosons. We will consider the density structure 
factor 

S{k,u) = J dxdt l ^ kx ^{p{x,t)p{Q,Q)), (21) 

and the spectral function A(k,w) = — ~lmG(k, w)signa; where the Green's function 
G(k,uj) is defined as [53] 

G{k, u) = -i [ dxdte*^*-^ (T$(x, t)^ (0, 0)]) . (22) 
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Figure 1. (a) Dynamic structure factor S(k, u>) and (b) spectral function A(k, ui). 
Shaded areas indicate the regions where they are non-vanishing. Liebs particle 
mode £i(fc) and hole excitation mode E2{k) are indicated. For spectral function 
A(k,cj) region with uj > 0(ui < 0) corresponds to the particle (hole) part of 
the spectrum. The figure above indicates the notation for the prefactor of the 
singularity in the response functions at a given edge of support. 



We consider the density structure factor S(k, u>) in more detail below, and present 
only final results for the spectral function A(k,ui). The exponents of S(k,u)), fj,i t 2 at 
Lieb's collective modes [38l [43] Si^(k) can be written as /Lt^a = 1 — fori — fi>L, where 
P'R(L) denote contributions from right (left) branches and are given by [29| . 

mL) ={— ± ^ + —) ■ (23) 

The phases 5± can be obtained explicitly by knowing the analytic form of the 
dispersion curve, the Luttinger parameter K and momentum k |29j . or by directly 
extracting them from the microscopic model (see Ref. [2TJ for the calculation of the 
phase shifts for bosons using the Bethe ansatz for the Lieb-Liniger model). 

In addition to the quasi-Fermi points, the field theoretical description of the 
singularities involves an impurity moving with velocity v d = dEx^fy/dk, and 
S(k, lj) in the vicinities of collective modes are written as 

S(k,u) =S h2 (k) J dxdte iSujt D(x,t)L(x,t)R(x,t), (24) 

where 5u> = lo — ei^fc), D(x,t) = S(x — v d t) is the impurity correlator, L(R)(x,t) = 
(i(vt ± x) + 0)~i 1l ( r ) , and we introduced prefactors Sx^k) which will be fixed from 
the comparison with form factors. In the vicinity of £2(k), Eq- (|24|) results in 

S { k,.) = 6(S.) , ^5 2 (fc)^-i 

v ' v 'r(fi R + fi L )(v + v d )^\v-v d \^^ ^ > 

while in the vicinity of £\(k) one has a two-sided singularity, 

. sm*PL L 0(5w)+Bmnji R O(-8w) 2nS 1 {k)5uj^+^ 1 

b(K,Lj) = ; ; — — ; r=— , (26) 

smvr^L + fin) T(jjk R + fx L )(v + v d )^ \v - v d \»* 
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In finite size systems, L(x,t) and R(x,t) get modified, see Eq. (fT2"]) . The change 
of D(x,t) to J2n D e 2 ™ nD ( x ~ Vdt ^ L corresponds to the quantization of the impurity 
momentum. Analysis of the scaling of the form factor with n r = ni = then leads to 

(27) 

where the exponents are specified in Eq. ([23"]). and \k;N) denotes a state of ./V bosons 
with a single particle (hole) with high momentum, and a hole (particle) at the right 
quasi-Fermi point, such that the total momentum is k. The state on the right is the 
full N particle ground state (in this case that of the bosons). We note that in Eq. (|2"T]) 
k has to be fixed before taking the limit L — > oo, since the k — > and L — > oo limits 
do not commute in the nonlinear Luttinger Liquid theory [23 [2H1 ■ It should also 
be noted that while our discussion focused on bosons in the Lieb-Liniger model, the 
relation in Eq. (|27|) is general and will apply equally well to the DSF of the CSM, 
albeit with different exponents fiR.L- 

Similar to the density structure factor, the spectral function also displays singular 
behavior near the Lieb modes [3TJ [55] and we find relations for the prefactors of 
the spectral function in terms of form factors of the creation/annihilation operators. 
We set the following notation before presenting equations similar to Eq. ([27]) for the 
prefactors of the spectral function. In the vicinities of £\(k) and —E2(k), respectively 
the spectral function behaves like 



with exponents [29] 

M± = 1 - Jm - W , 

— _(Vk s ± (k)Y 

- [~2 — ^r) • (29) 

Note that in the above the phase shifts 5± can be calculated by knowing the dispersion 
curve and k [55] or from microscopies [2T] . In Ref. [3T] , the phase shifts are bosonic and 
calculated from the Bethe ansatz for the Lieb-Liniger model. In Eq. ([2"9")l above, the 
expressions for the exponents refers to fermionic 5±, consistent with the notations of 
Refs. [21] [15]. The final answer for the exponents in Refs. [21] [29] are all in agreement 
provided the correct 8± are used. 

By analyzing the scaling of the form factor of the creation/annihilation operator 
with n r = m = we obtain 

|< t ,„ Tl |^| W) |^«(|p\ (30 ) 

with /TrTJ) specified in Eq. ([29]) . The state \k, JV±1) refers to a state with an additional 
particle or hole such that the total momentum is k, while \N) is the ground state of 
N bosons. 

Similarly, near S2(k) and — £\(k) we have 

T(^i R + fj, L )(v + v d )^-\v - v d \i^- 
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with exponents [2"9"] 

H± = 1 - jJ,R - flL, 

(SVK 1 6 ± (k)\ 2 

We note again that the phase shifts 6± are meant to be fermionic in line with notations 
of Refs. [351 US]- Alternate expressions for the bosonic exponents with bosonic phase 
shifts calculated in terms of microsocopics of the Lieb-Liniger model can be found in 
Ref. [21]. Moreover for the prefactors we obtain 

|(fc,7V T l|V< (t) |A0| 2 «^^ T ] , (33) 

with /ifl(L) specified in Eq. (l29l) . The state \k, N±l) refers to a state with two particles 
(holes) near the right quasi-Fermi point and an additional hole (particle) such that the 
total momentum is k, while \N) is the ground state of N bosons. Let us also mention 
that for certain parameters A(k,w) can have a non-analyticity at —e\(k) instead of 
a divergence [29]. In this case A_(k) refers to the prefactor of the non-analytic part. 
The results for prefactors of the singularities in the spectral function for fermionic 
models are presented in Ref. |T5| along with expressions for the relevant exponents. 

Finally let us note that the methods described here can be also used to fix the 
prefactors of dynamic response functions of integrable lattice models, such as the 
model of spinless fermions with nearest neighbor interactions [25] or XXZ model. For 
that, one needs to combine their form factors [48] [70] with the corresponding field 
theoretical description of the singularities in their DRFs [24[ [25l |34] . We postpone the 
discussion of DRFs of these lattice models for future due to more complicated nature 
of their exact eigenstates. 

3. Density structure factor of the Calogero-Sutherland Model 

We demonstrate here how the relationship between non-universal prefactors of 
correlation functions and form factors of excited states calculated for the corresponding 
operators, established in Eqs. (111)1) and ([2"?]) can be used to calculate prefactors in the 
correlation functions of the Calogero-Sutherland Model (CSM)[3S1 [3B1 [3Z] described 
by the Hamiltonian 

ft - V 02 + V MA - 1)*' 

Here N is the total number of particles of mass 1 /2, contained in a system of length L 
held at density p = N/L and A = p/q > 1/2 with p, q coprime controls the strength 
of the long range interactions between the constituent particles. 

The CSM admits an exact solution and consequently one obtains the energies 
of the ground state and the spectrum of excitations. Furthermore, the CSM ground 
state wavefunction has a product form |35) 

^GS = Y[(Zr ~ ^) A rK A(JV ~ 1)/2 ' Where Z l = ( 35 ) 

i<j k 

and can be interpreted qualitatively as the wavefunction of a gas of non-interacting 
particles with fractional exchange statistics. The excited states can moreover be 
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constructed from the ground state by multiplying by Jack symmetric polynomials |36j 
and are completely characterized by a set of N quantum numbers n\, ...,n^. These 
quantum numbers allow us to solve for the "asymptotic" quasimomenta which allow 
us to calculate physical observables like energy and excitation spectra, using 37 : 

27m,- tt(A- l)(2j - N - 1) , . 

Pi = -T + - J = h-,N. (36) 

One may use the intuitive appeal of asymptotic plane wave states to interpret the 
quantum numbers rij as the wave numbers of these asymptotic states. The algebraic 
structure of these polynomials is efficiently captured in terms of operations performed 
over Young tableaux [7T] [75] . 

The CSM exhibits an additional property that leads to the existence of closed 
form expressions for form factors of physically relevant operators. The action 
of, e.g. the density operator, on the ground state can only result in particular 
types of states allowed under "selection rules" , corresponding to states with a finite 
number of excitations [45]. Such a structure is a consequence of the fact that 
the ground state wavefunction and the operators acting on it admit symmetrized 
representations in terms of Jack polynomials which have to satisfy orthogonality 
relations. Qualitatively, creation of states with finite numbers of excitations can be 
thought of as a generalization of the action of the operator on the gas of free fermions 
where it results in a finite number of (real) particle-hole processes. Thus, the space 
of eigenstates of the CSM possesses a structure reminiscent of Fock space, which can 
equivalently be seen from the form of the ground state wavefunction. This property 
greatly constrains the complexity of the form factors of the CSM and allows one to 
obtain analytic expressions. 

An expression for the density-density correlator is available for a finite sized 
system for the CSM due to Ref. [IS] . One begins with a representation of the density 
fluctuation operator as 



1 N N 1 / °° N 

2=1 \m=l 2=1 



e-^'t + c.c. , (37) 



which may in turn be expanded in terms of Jack polynomials. Using the orthogonality 
relation for these polynomials the density-density correlator may be written as a sum 
over Young Tableaux [45]: 

<P(*,*)P(0,0)) = ^E ^^^^ (38) 
where in the above expression the following notation is employed: 



/r, \2 N 

H" = II (a+(i-l)A-(i-l)), 

(-8,2) ere 

A - II U-, + ^p)U-, + i + ^ 

(i,2)ere V 7 V 



(39) 



Exact prefactors of correlation functions of ID quantum integrable models 



11 



where Ki and n'j refer to the total length of the i th row and the j th column respectively. 
Moreover |«|, called the "weight" of a partition is the total number of blocks appearing 
in a given Young diagram. The Young diagrams above are indexed by k and particular 
"cells" or blocks contained in the diagram are indexed by (i,j), with i,j > 1. Note 
that the prime in [0'] is an instruction to skip over the cell i = l,j = 1 which will 
cause all terms to evaluate to 0, see Eq. (|39|) in the product. 

The expression in Eq. ()38[) is in direct correspondence to the form factor expansion 
of the correlation functions, see Eq. (fTU|) . which involves a sum over all excited 
eigenstates of the system; every term in the sum Eq. (f3"5|) explicitly contains the 
momentum k and energy w dependence of the state in the exponential term. Thus we 
may uniquely identify the energy and momentum of the state associated with a Young 
diagram in the sum. This provides a recipe for constructing Young diagrams associated 
with excited states created by the action of the density operator on the ground state; 
we start by ordering the quantum numbers of the ground state and record the shift 
in the i th quantum number of the excited state as the length of the i th row of the 
Young diagram («j in the above notation). It should be noted that the selection rule 
mentioned earlier is encoded in the term [0']^ which appears as a coefficient when 
expanding the density operator in terms of Jack polynomials in Eq. (|37j) . This term 
is zero unless the Young diagram n contains a p x q block. We may interpret this 
as an excitation of p quasiparticles and q quasiholes, and the mathematical structure 
reflects the fundamental way in which the density operator may act on the ground 
state. 

We may directly obtain the form factor contribution from an excited state with a 
definite energy and momentum in terms of algebraic operations on a Young diagram. 
However to extract the prefactor of the correlation function in the thermodynamic 
limit from the form factor, one needs to carefully separate the non- trivial power-law 
of the system size L, see Eq. (|27|) . from the remaining contribution. We carry out this 
analysis below for the prefactor of the density structure factor and obtain analytic 
results. 

We wish to identify a single diagram that captures the contribution to the density 
structure factor when lj « 1^(^)1, i-e. in the vicinity of the edge of support with 
e{k) = X(k 2 — hp), and evaluate the corresponding form factor in the thermodynamic 
limit. We know from the finite size analysis of the effective three subband model, see 
Refs. [HI [29], that the contribution to the prefactor iS'(fc) comes from the form factor 
of the lowest energy state of momentum < k < 2kF- Such a state is the CSM analog 
of free fermions with a particle near the right branch, and a hole of momentum kp — k. 
Following our prescription for constructing Young diagrams for states with particle- 
hole excitations over the ground state, this corresponds to a Young diagram with a 
single column of length c\ = ^-k — (p — l)(g l — 1). We note here that there are also 
singularities for k outside the range (0,2/s.p) and the prefactors of these singularities 
can be obtained from form factors of the state with additional umklapp excitations (a 
hole on the left branch, a particle on the right) on top of the configuration described 
above. For simplicity we present only the result for the first prefactor. We will 
choose groups of terms appearing in the form factor expansion in Eq. (|38j) to evaluate 
together. This is necessitated by the fact that some terms may individually diverge 
factorially (~ N N ) in the thermodynamic limit, but their divergence is suppressed 
by other terms in the expression in such a way that the full expression will contain a 
non-trivial power law in L and a prefactor. 
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(a) 



basic block 
V 



(b) 



T2 



1 n 



quasiparticles 




quasiholcs 



0-2 



Cl 



Cl = ^fc-(p-l)( g -l) 



Figure 2. (a) Typical Young diagram which appears in the sum in the expression 
for the density-density correlator see Eq. 13811 . The basic block appears in all 
non-vanishing contributions to the density-density correlation function and is a 
q X p block. Excited states correspond to states with quasiparticle (quasihole) 
excitations, rows (columns) of total length ri,r^, ...(ci, ca, ...) above the basic 
configuration, where limjv— >oo c,r/N ^ 0; (b) The Young diagram corresponding 
to an excited state of momentum k, with energy e(k) = A(fc 2 — k 2 F ) which defines 
the edge of support. For free fermions this state would correspond to a state with 
a particle-hole pair with momentum transfer k on top of the N particle ground 
state 



We first consider the following group of terms: 
r _ 



(40) 



Let us first evaluate the contribution from the basic block (see Fig. ([2])) which 
should be present in all terms: 



q p (n+ ^ - (i - 1)) 

rr = niP 



=1 j=l 

1 



N- 



n 



XN + j-1 



n 



N + l-i 



N J 1L yXN-p + X + jJiiyN + q-i 

j=l v 7 1=1 v ' 

= 1 + 0(1/N). (41) 
Additionally, the presence of particle-hole like excitations i.e. columns with length 
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large in the thermodynamic sense, will lead to the following contribution: 

N-i + 1 



tt= n 



l=q+l A 

r(JV- g + i)r(JV-ci + i/A) _ / k 

T(N-q + 1/X)T(N - ci + 1) ~ V 2fc F 



(42) 



Next we consider the term T2 = ■ 



We again calculate separately the contribution from the first column 
corresponding to the particle-hole like excitation: 

r A (*-d 2 ft (*-d 2 

^ 1 + 1^ 1 + 1 + « - VA) (ci - i + 1/A)(d 

i\, r( Cl ) 3 r( Cl + 1 - i/A) 



A y r(ci + g)r(ci + g + 1 - 1/A)r(ci - g + 1/A)r(ci - g + 1) 



2tt J \\ 

The remaining terms evaluate to unity when properly normalized by requiring the 
basic block form factor evaluates to po- We may now combine the above expressions 
to obtain a prefactor of the DSF using Eq. (|27J). the correspondence between the 
Luttinger parameter K and the interaction parameter A, j = K, and the exponents 
Ail, P-L given by 

^ , , , 1 , 5 +( _ } (fc) n2 



1 \ 




2V^f 



(44) 







t L 


lim 


2ir 




L—too 




K 2^ 





















|(A;;7V|p|0,7V)| 2 



2fc F fc 



2fc, 



(45) 



We were also able to evaluate the equal-time density-density correlation prefactor 
A\, see Eq. ([1]), by considering the form factor for a state with a single column of 
length TV over the basic configuration. We obtain for the prefactor 

A, = 2 r(l + ? )2 , (46) 
(2tt) a 

in agreement with the result of Ref . [12] , where general expressions for the prefactors 
of oscillating components of the density-density correlators were obtained using the 
Replica Method. 
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4. Prefactors of the Lieb-Liniger Bose Gas 

The one dimensional Bose gas with contact interactions of strength c > is described 
by the Hamiltonian 

H = J dx [d x ^ {x)d x ijj{x) + c^{x)^(x)^(x)^p(x)] , (47) 

where ip',i/j are the boson creation and annihilation operators, and where we assume 
the particles to have mass 1/2. The problem has an exact solution given by the Bcthc 
Ansatz [38] . The solution of the A-particle interacting Bose gas is obtained in terms 
of well defined ground and excited states characterized by sets of N quasimomenta 
{A} and {/i}, respectively. 

The quasimomenta of the ground state for the problem are given as the solutions 
of the Bethe equations [33] : 

N 

L\ 3 + - X k) = 2 ™j = M3 ~(N+ l)/2), j = 1, 2, N, (48) 

where 9(X) = i log ^ffrj^ = 2 arctan (~) , and where the equation carries the physical 
meaning of a particle traversing an entire turn of a ring and returning to its origin 
having picked up a total phase shift resulting from pairwise two particle scattering 
from the other particles in the way. 

An arbitrary excited state may be created from the ground state by simply erasing 
some finite to < A of the nj, {n^[ , , n^}, and replacing them with new quantum 
numbers {nf , n% , n+ }, that are not from the original set of rij. If we call the 
new set of quantum numbers {fij}, then {fij} is obtained as {nj}/{n~} U The 
quasimomenta fij characterizing this excited state are given by the solutions of the 
equations: 

N 

Ln J + Y^(Pj-M)=2imj, j = 1,2,..., A. (49) 
fc=i 



4--1- Form Factors 

As demonstrated in the previous sections, full knowledge of matrix elements calculated 
from the microscopic theory is sufficient to obtain the associated prefactor of a 
correlation function. Thus, we wish to obtain the thermodynamic limits of various 
form factors between ground state and excited states specified by the field theory. 
The final expression for the limit of the form factors is expected to have a power law 
behavior as a function of the system length L, see Eqs. ([T5 |) - (|16p . and Eq. (|2"T|) . Again 
our task is to separate the power law from the prefactor in the form factor when faced 
with terms that are badly divergent (0(N N )), so we proceed in a similar spirit to the 
calculation of the DSF prefactor of the CSM. 

We begin by considering the exact expressions for the form factors of the finite 
size Bose gas that arise from the machinery of the Algebraic Bethe Ansatz [43] . 

The form factor of the density operator between normalized eigenstates |{/^}jv) 
and | {A} at) is given by the following expression [46] : 

"IPex 



^=({/i}iv|p(0)|{A} JV ) 



exp[-ixP ex ] - 1 



(Mat I / dyp(y) | {A}^ 
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where 



Xjk + ic\ Det(Sjk + U 3 k) 



(50) 



n 



Hk - Aj ± ic 
At. — A,- ± zc ' 



(51) 
(52) 



Ujk - v+ v- W 



Am A i 



X(A j -A fe ) = J ftT(A i! A fe ) = 



(^(Aj-A^-^Ap-Afe)), (53) 
2c 



(A,- - A fe ) 2 + c 2 



(54) 



Note that in the above expression p is any integer and this freedom of choice is due 
to the rank deficient structure of Ujk 06], [65] . 

Similarly, the form factor for the boson annihilation operator is given by [46] , [65] 



\({v}n-i\iP(0)\{X}n)\ 



IK/4tf-i|| 1/a l|{A}*ll 1/2 



(55) 



where 



JV-l 

1 j,k j 



.1/2 



N-l 



n & 

j,k,k^N 



-Det JV _ 1 (5 jfc + %), (56) 



JJ Hk - Aj ± ic 



fc=i 



A' 



JJ Afc - A. ± ic 



fe=i 



-A,) 



v- + -^n^(A m -A,-) 

Moreover, the form factor for the creation operator Q + is obtained from the one 
for the annihilation operator Q~ by Hcrmitian conjugation. 

In the above expressions the norms are given by [58], [73], [46], [43] 



(57) 



(K(Xj - A fe ) - K(X N - A fe )). (58) 



!{A} A : : =c^JJ Ai " 



Xj — Afe 



Dct 



N I $jk 



N 



- K{Xj - Afe) 



(59) 

Using the density operator and boson creation and annihilation operators, we may 
obtain expressions for the density structure factor (DSF) and spectral function (SF) 
using the Lehmann representation [66) . 
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4-2. Thermodynamic Limits and Finte Size corrections 

The various terms in the expressions for the form factors, Eqs. (|50)l - (f55"j). depend 
directly on the quasimomenta of the ground and excited states. A thermodynamic 
description should not make reference to these discrete quasimomenta. We wish to 
define "thermodynamic quantities" that provide an equivalent, but sensibly defined 
description in the limit of large N,L. In the thermodynamic limit, the ground state 
of the Bose gas can be described using a continuous distribution of quasimomenta A 
belonging to the set of real numbers (—q,q), while some excited state with n holes 
and m particles consists of quasimomenta in the set (— q, q)/(X^ . ...A~) U (A^, A+). 
For this purpose we will use the density function p(X) and the shift function F(A|^), 
that are conventionally used in thermodynamic descriptions of the Bose gas. These 
functions are given by [43] 

p{X) ~h[_ W A ^)P(M) = (60) 

nm - ^ £ dnK{\ - p)foi\v) = (6i) 

where the shift function picks up an overall positive or negative sign depending on 
whether the excitation v in Eq. ([BTj) is a particle or hole. Thus, the shift function for 
a particle-hole pair excitation can be described by F(A|/x + , A - ) = F(X\p + ) — F(A|A~). 
Moreover using the linearity of Eq. ([BTj) . a general excited state with multiple 
excitations of both types can be written as a sum of shift functions taken with 
appropriate signs. 

Furthermore, we may express the number density, poi 01 the Bose gas using the 
ground state density function 

* dXp(X) = po. (62) 
'-<? 

To successfully account for all constant factors in the thermodynamic limit, we 
will need to consider finite size corrections (in orders of 1/L) to these functions and 
collect any corrections that sum to finite values in the thermodynamic limit. We 
begin by writing the system of N Bethe Ansatz equations that give the ground state 
quasimomenta of the N particle interacting Bose gas 

N 

LXj + £ ~ A fe ) = Tmrij = 2n(j - (N + l)/2), j - 1, 2, N. (63) 
fe=i 

The quantum numbers rij on the right hand side give us a natural way to index the 
quasimomenta as we take the above system over to the thermodynamic limit - we 
will want to keep the ratio of rij/L a constant as we take L — > 00. It is convenient 
to construct a generalized continuum version of the n^s that smoothly interpolates 
between them in the thermodynamic limit, as described in Ref. [43 . 
Let us define a variable x that satisfies 

N 

LX(x) + - Afc) = 2Lnx. (64) 

k=l 

We now define a function pl as 

Pl(X(x)) = Jg. (65) 
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This function carries the usual meaning of the density of quasiparticles in a given 
quasimomentum interval, but we will require it to retain finite size corrections. 

To obtain finite size corrections to our thermodynamic quantities we may use the 
Euler-McLaurin formula |74) which quantifies the difference between a sum and an 
integral in terms of a series in powers of the discretization (in our case 1/L). We 
obtain for pl [43] : 



1 

~ 2^ + ~ 
1 

~ 2^ 
1 

~ 2tt 



1 

2^ 



1 



N 



K(X(x)-X k ) 



k=l 



1 



2ttL 
1 

2^ 
1 

2^ 



N/2L 



K(X(x) - p(y))Ldy 



-N/2L 



1 B 2 
2ttL 2 



^dK(X(x 



Mi/)) 



Ldy 



V—N/2L 



y=-N/2L 



K(X - p)p L (p)dp 



1 1 fK'(X-p) 



K{X - p)p L (p)dp - 



2ttL 2 24 
1 



(K'(X-q)-K'(X + q)) + 0(1/ L 3 ). 



48nL*p L (q) 

Let us note two things about the above calculation. Firstly, in the above 
expression the discrepancy between the boundary value of x b . c . = N/2L and nj^/L, 
the quantum number associated with the edge of the distribution q can be resolved 
by modifying the definition of q to absorb the linear term from the Euler-McLaurin 
expansion [43], i.e. we redefine q as, q — > Ajv + 2Lp(\ N ) • Secondly, we will only need to 
keep track of corrections to 0(1/ L) in the subsequent calculations since this captures 
all finite results in the thermodynamic limit. 

We would like to characterize the excited state with the same precision as the 
ground state. Excited states are obtained from the ground state by replacing a finite 
number of ground state quantum numbers rij with new quantum numbers hj that are 
not from the set of ground state rij . For a choice of finite N, L these excitations are well 
defined in that we will have two finite sets of quasimomenta {p + }, {A~} corresponding 
to the particles and holes respectively. The notation is intended to be suggestive of 
the fact that p + are "actual" quasimomenta in the excited state, while A~ will occur 
only as a subset of the set of ground state quasimomenta and are absent in the excited 
state. 

It will be more convenient for our purposes to also consider "artificial" hole 
quasimomenta {p~} and attribute them to the excited state. This is meant in the 
following sense - all the quasimomenta of the excited state are shifted with respect to 
the ground state quasimomenta. We expect this shift to be of 0(1/ L). Had the hole 
quasimomenta been present in the excited state, they too would have been shifted with 
respect to their ground state counterparts, {A~}, on the order of 1/L. Consequently 
we may realize p~ as A" + 0(1/ L). 

We would now like to obtain thermodynamic versions of {p },{A - } and retain 
0(1/ L) corrections to them. 

Let us start by considering the simplest case of a single particle-hole pair and 
some finite N, L. Furthermore let us denote the quantum number corresponding to 
the particle, n + , and the quasimomentum, p + , with the corresponding hole and its 
counterpart in the excited state defined as A - , p~ respectively. 

For convenience, let us define the following notation. A single prime 
accompanying a sum, , means we are leaving out terms corresponding to p + , A~ in 
the sum. A double prime accompanying the sum, ^ , means we are adding in p~ in 



(66) 
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place of p + . To go from ^ to , we explicitly treat the finite number of excitation 
terms. To go from to J^" we compensate the extra terms with terms outside the 
sum of the opposite sign. 

Starting from the Bethe equation for the particle excitation we proceed as follows: 

N 

Lfj+ = 2im + - #0 + - Mfe) 

k=l 

Lp+ = 2vm + - ' 6(ji + - p, k ) - 6(p + - pT) + 6{p, + - pT) 

N 

Lp+ = 2vrn+ - ^((?( M + - X k ) - K{p+ - X k )(p k - A*)) + - A") 
fe=l 

-i<r(M + - A~)0*- - A") + 0(1/L 2 ) 



A' 



'0Q i +-\ k ) + K(i*+-\ h )%$\ 9(fi+-\- 



m + = 2 -ir - E I — — — — f + ^ L ^ ; + o(i/l 2 ) 

p+ = 2n— - dX9(p + -X)p L (X) + -l- dXK(p+ ~ X)F{X)+8(p+ - X~) 

+0(1/L 2 ). (67) 

Note that in the above derivation we have used the fact that we expect p k — X k ,p~—X~ 
to be 0(1/L), and the relation, (Xj — pj) = Lp(\\ ^(V^ 2 )- Using the last step 
of the above as our starting point, we may now take the thermodynamic limit. The 
prescription we use is to keep n + /L constant as we send N, L — > oo. Furthermore 
let us separate the 0(1) and 0(1/ L) contributions to p + for the sake of clarity, i.e. 
p + = p% + [xf, L + 0(1/L 2 ). We obtain 

i4 + i4/l = 27r — - / dxe ^i - x )p( x ) - i4/l j dXK (ri - x )pw 

q dXK(p+ - X)F(X) + 9(/4 - X-)) + 0(1/L 2 ), 

(68) 

where on the right hand side we have substituted p for etc. because the lowest 
order difference between such terms is higher order in 1/L than we are keeping. 
Thus we obtain as the thermodynamic limit and the first order correction 

M + = 27T— - / dX9(p+ - X)p(X), 



L 



-q 



+ _ 1 / Jl q dXK(^+ - X)F(X) - e(^+ - x-) \ _ 
th/L L { 1 + j q _ q dXK{p+ - A)p(A) ) Lp(p+) ■ 

(69) 

In the last step above, we have used Eqs. (|60p. (|6"Tj) . 

Similarly, let us consider the Bethe equation for the ground state quasimomentum, 

x-, 

N 

LX- = 2nn- - ^ 8{X~ - X k ) 
fc=i 
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„- ri 

\-=2ir—- dXO(X~ - X)p(X) + 0(l/L 2 ). (70) 

L J-q 

We wish to now obtain pT . Since we expect p~ = A - + 0(1/ L), and the lowest 
order correction comes from the usual shift in the excited state quasimomenta, we can 
quantify this exactly as 

ri 



nT f q 
Po = 27r — - J d\9(n - X)p(X), 



Ml/L ~ Lpfay (71) 

We are now ready to obtain the lowest order finite size corrections to the shift 
function. Let us start with a definition for the finite sized shift function. For an 
excited state with some set of particles and holes, {/i^,/!^, ...,p~}, 

where pT is used in the sense described above, we have 

Fl{W{i4,i4 >-»Mn}W.M2 >-.Mn» = (Aj - Pj)Lp L (Xj), (72) 

where Aj , pj are the ground and excited state quasimomenta respectively. 

We may obtain an equation for Fl by subtracting the Bethe equations for pj and 
Aj for jjbj not in {p + ,p~}: 

N 

L(\j - Pj) + J>(A,- - A fe ) - 9{jij - p k )) = 
fe=i 

L(Aj - /ij) + ^2'(0(\j - A fe ) - (9(jUj - p k )) + 0(/ij - p~) - 6{nj - p~) - 9(pj - p + ) = 
l ( x j - Pj) +^2 {6(Xj - Afe) - 6(Xj - Afe) - K(Xj,X k )(nj - p k - Aj + X k ) 
-\ K \ X 3 - Afe)(/"j - Mfc - Aj + Afe) 2 ) = 9{nj - /i + ) - 0(/Zj - p~) 
L(Xj - (l + j Y,"K{Xj,Xk)) = ^2" K(Xj, Xh)(X k - p k ) + 9{X j - f4) - 9(Xj - 



+ ^2"^K'(Xj - X k )(nj -Hk- Aj + Afe) 2 +K(Xj,i4)(nj - Aj - i4/ L ) 
-K{Xj,^){nj -Xj-fj,- /L ), 

2nL(X J - pj)p L (Xj) = 52"K(\i, X k )^-^-Lp L {X k ) + 9(Xj - p+) - 9{X j - ^) 

LpL\Xk) 

+ ^2"^K'(Xj - X k )(nj -Hk- Aj + Afe) 2 + K(X j ,p+)(p J - Aj -m+ /l ) 
~ K (Xj,Po ){Pj - Aj -Mvl) 

*i(Mj/)) 



/IV / 
LdyK(X,p(y)) 
-N/2L 



L£^ ( , mA _ rt (^_^l) 2 +0(1/l2) . 
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1 



2ttL 
1 



2n 



+ , (F L {\) F L (p+) 



-^(A,Mo) 



fF L (X) Fl (pp~ ) 
\Pi( A ) Pl(Pq) 



dpp L (p)K'(X - p) 



Fl(X) F L (p) 
Pi (A) pl (a) 



0(1/X 2 ). 



(74) 



Note that in the above derivation, Fl(X) is a condensed notation for Fl(X\p + , p ), 
and moreover this Fl is related to the more general one defined in Eq. (|72|) as 

Fi(A|{/i+ /i+ ..,/i+};{/ir 5 M2>-: Pn» =F L (\\p+,p-) + ... + F L (\\p+,p-). (75) 
Furthermore, we expect that in the thermodynamic limit Fl(\\pI , p^) = 
F(A|^+) - F(X\pr) + 0(1/L), where F(A|/if ) is given by Eq. (EE]). 

It is also convenient to define an analog of the ground state momentum density 
function for the excited state. This will help in calculating the normalization of the 
excited state that appears in Eq. ([5U|) and Eq. ([55| . Before proceeding, let us write 
down the integral equations for the derivatives of pi, , Fl as these will be used in the 
simplifications to follow. 

ri 



2nF' L (X) 



dpK'(X - p)F L (p) + K{X j - p+) - K(Xj 



Po 



-9 
'/ 



2tt Pl (X) = \ dpK'(X - p)p L {p) + 0(1/ L 2 ). 



0(1/L), 
(76) 



To obtain the density function for the excited state, we start with the Bethe 
equation for a generic quasimomentum from the excited state, Eq. (|49[) . and 
differentiate once with respect to pj : 

L = 2irLp ex , L (pj) - ^ K iN ~ Pfe) _ K iN ~ M~) + K (l i j ~ lO _ k (Pj ~ P + ) 

k 

L = 2TTL PeX:L (pj) - ^2 K (Pj - Pfe) - K ip-j - P + ) + K{H - P~) 

( F L (X 3 ) F L (X k ) 



k 
N 



L = 2irLp eXtL (Xj) - 2J 



fe=i 



K(X 3 - A fe ) - K'(Xj - A fc ) 



-2VL(A,)T#|-iC(A, 



Pl{Xj) 



2np ex (X)-2np' L (X) 



Fl(X) 



\Lp L (Xj) 
ti)+Ki\ j -to) + 0(l/L) 
_ 1 K(X-p+)-K(X-po) 



Lph{X k 



dv 



L P l{X) 

K(X-v)p L (v)-K'(X-u) [ p L (v) 



L 

Fl(X) 
L P l(X) 



Fl{v) 



0(1/ L 2 ) 



2tt Pl (X) + 2ir6p{X) - 2tt P ' l (X) f L \}. = 1 + / dvK(X - v)p L (v) 



Lp L (X) 



-2n P ' L (X) 



Fl(X) 



-2tt 



F' L (X) 



(77) 



L PL (X) ' " L 

The first two steps of the above set of equations is exact. From the third step 
on we retain only terms up to 0(1/ L). In going to the third step we have expanded 
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all terms dependent on p,j in terms of Xj as this is what is useful for us in our final 
expression (see the following sections). In going to the last step we have used the 
relations given by Eq. ((75)) . 
From Eq. (|77|) we have 

p ex M - Pl{X 3 ) = ~ (fl(Aj) - ^jffi ) + °(V^ 2 )- (78) 

It is useful to consider the notion of a "partial thermodynamic limit" . For instance 
the algebraic expression for the norm of the eigenstate appearing in Eq. (|59|) can be 
expressed in terms of a partial limit as established in Refs. [43], [58], [46] 

lim ||{A}|| 2 = J] ^±^Y[(2nLcp(X,))Bet (l - ±-k) , (79) 

3 >k 3 k j v 7 

where the determinant is intended to be a Fredholm Determinant (see Ref. [75 , 
Appendix D for more details). The above expression is a "partial" limit in the sense 
that it contains both the quasimomenta of the ground state which should be absent 
when the "full" thermodynamic limit is taken, as well as the Fredholm Determinant 
which is well defined in the TDL. 

Similarly, we may write down the norm for the excited state as 

IIM " 2 = n ^ n^^OPet (l - ±k) . (80) 

j>k 3 K j x ' 

We will frequently use the notion of such partial limits because individual terms in 
the algebraic expressions for the form factors may not have good thermodynamic limits 
- they may often diverge in a non polynomial way. However, they can be regrouped to 
obtain expressions that have good (finite or power law divergent) limits. Thus we want 
the final expression for the form factor to be well defined in the limit, on the other 
hand the intermediate steps leading up to this will be a mixture of thermodynamically 
well defined quantities and the aforementioned partial thermodynamic limits. We will 
take partial limits of simple groupings occurring in the form factors in the following 
sections and later collect them and express the final answer for the form factors. 

Let us consider the term /V~ with Vp as defined in Eq. (|52[) 
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Now using the relations given by Eq. (|76[) we have 
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Substituting in Eq. (|81[) we obtain 



(82) 



log 



Yl 

v; 



2mF L {Xj) 



1 



1 



Lp L {Xj) 



F' L {X 3 ) 



l l.(Xj V'/ / A ; 

2px(Ai) 



0(1/L 2 ). (83) 



Thus we may handle the following term appearing in 57, Eq. (|5ip . as follows: 
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where in the above expression the term <^>(j) is purely a phase. Such terms will drop 
out in the final answer because we only need the absolute square of the form factor. 
The product of sine terms can be expanded in the following way 
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(85) 



The other terms in fi may also be regrouped 
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With this form for fi we may combine it with the norm terms : 
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(87) 



In going to the last line we have collected a purely phase term in the last product. 
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After substituting the expression for V and regrouping terms we obtain the final 

form 



|F| = 



1/2 

-q i (Xj - p k + ic)(\j - /<tfc + ic) x 



n 

i=l 



(Xjk + ic)(fijk + ic) 

F L (Xr) \ iP ex Det(5 3k + U 3k ) 



1 /2 

tt sin(7rF L (Aj)) -pr / X jk [i jk \ 
11 ttF l (A,) lllfe-A fc ) 2 



L( PL (A 4 )p ea; (M+)) 1/2 (A,; - (7+ - 7 p -)Det (l - f ) 

^F i (A)cos(^F i (A)) ^ (A) _ F L (AK(A) \ _ 1 / , (A) _ F L (A)p' L (A) 



° M> ' ' sin(^F L (A)) V^" V ' V 2 PL (A) ; 2 1^"' Pl(A) 

Similar manipulations can be performed to obtain the partial thermodynamic 
limit of the form factors associated with the boson creation/annihilation operators. 
For the creation/annihilation operators we consider a matrix element between states 
with different numbers of particles (see Eq. (]5"5"]) ). For the creation operator form 
factor, the ground state has N particles while the excited state has N + 1 particles, 
while for the annihilation operator form factor the ground state has N particles while 
the excited state has TV— 1 particles. This means that in the construction of the excited 
state, there is either a missing or extra quantum number, relative to the ground state. 
This difference should manifest as an extra hole (for annihilation operator) or particle 
(for creation operator) in the excited state. To account for this key difference from 
the calculation of the density operator matrix element, we will introduce modified 
shift functions to express terms in the form factors of the creation and annihilation 
operators. We define these modified shift functions as 

■ + ,~,l 

,+ 



F_(A) = F(A|K, {hi, m" +1 }) - tt P (A), (90) 

for n > and with F + appearing in the creation operator form factor and F_ appearing 
in the annihilation operator. The terms in Eqs. (|89[) . (J90J) can be justified as follows. 

The first term in Eqs. (|89|) . (|90| corresponds to the particle/hole excitations of 
the state including the extra particle or hole excitation due to the difference 
in number of particles in \fi) versus |A). This term is analogous to the shift 
function used for the density operator in Eq. ([75]). F(X\{fif, {fJ-i , M^}) = 
F(X\fil) — F(X\fi± ) + ... + F(A|^+) — F(A|^i~) where [if and \i~ are particle and 
hole quasimomenta with respect to the N particle ground state. The second term, an 
additional factor of ±irp(X) has a more subtle origin - the derivation of the equation for 
the shift function, Eq. (|61[) from |43j , assumes a change of periodic boundary conditions 
to antiperiodic ones with the change of particle number. Phase shifts calculated from 
this shift function correspond to the fermionic Cheon-Shigehara model [76] dual to 
the bosonic Lieb-Liniger model. To obtain the desired phase shifts, one may "correct" 
the shift function by adding Ttp(X) to the unpaired shift function in Eq. (f89| and 
Eq. (|90|) [21] . The sign before this correction is picked up according to whether the 
unpaired shift comes from an extra particle or hole. The reason these two terms are 
simply summed to give the correct total shift functions comes from the fact that we 
can sum the integral equation for the shift function for two different particle-hole 
configurations and obtain the integral equation for the composite shift function, i.e. 
the shift when all the individual particle hole configurations occur together. 

Furthermore, there are a few modified terms in the form factor of the 
creation/ annihilation operators that need to be carefully accounted for. To do so we 



.(88) 
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will use the following procedure, focusing on the annihilation operator for concreteness. 
Let us choose one of the n holes in the excited state and consider its quantum number. 
There should be a ground state quasimomentum corresponding to this quantum 
number, with some index which we will call n~ . We are assured of at least one such 
hole since there is always one less quasiparticle in the excited state used for obtaining 
the creation operator form factor. We will seemingly treat this hole specially, but 
this is just a convenience as we could have picked any of the potentially n available 
holes - but we will commit to one to simplify calculations. The final answers for the 
prefactors will be independent of the particular choice of n~ . 

Now we reindex the excited state quasimomenta as follows. For i < n~ — > fii, 
for i > n~ , fj>i — > fu+i- Thus, in all subsequent steps there is a "gap" in the indices 



F -(K-) 



Thus we have n 



1 



for fj, at )Jb n -. We will artificially define p n+1 = A n - Tp(\ — ) 

holes and n particles in the excited state with quasimomenta, /if for n > and 
we will "remember" the index of one of the holes, n~ . Note that fii is simply a 
reindexing of the original set of fii with no change in numerical values implied. Under 
this reordering it becomes possible to define F_(Aj) = (Aj — fii)Lp(Xi),i ^ n~ and 
F_(X n -) = (A„- — p~ +1 )Lp(X n -), as in Eq. (|72")l . where the F_(X) is the one defined 
in Eq. 

Let us consider the combination V + jV~ appearing in the form factors in Eq. (|57p 




F-AXM^y 

2pl(A,) 



0(1/^). 



Lp L (Xj) 

We use the double prime notation as described in earlier sections (see the text 
prior to Eq. (f67|) ) and collect the hole created by having one less particle in \n) with 
the other n holes of the excited state and raise the index for such terms to n + 1. To 
obtain the last line we have used the same steps as in Eq. (fSTj) . 

We will now consider the following grouping of terms appearing in Eq. (|55p 
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Moreover, the following terms from Eq. (|55p give 
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On combining these results with the rest of the annihilation form factor we obtain 
as the partial limit of Eq. ([55} 
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The derivation of the result for Q + is almost the same as the one for Q , hence we 
present only the answer: 
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Note that in Eq. (|94|) and Eq. (|95|) the notation Sjj, is used to differentiate between 
the matrices appearing in Eq. (|55|) when the excited state has N ± 1 quasimomenta. 
The matrices are given by 
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Lm^jV^m - Aj) 

It is more convenient to separately evaluate the thermodynamic limits of specific 
groups of terms appearing in the expressions Eqs. (j8"8"j) . (|M|). (|93|) . We begin with the 
terms in the square brackets. 



4-3. Evaluation of Mi 

We will denote by Mi the first group of terms appearing in square brackets in 
Eqs. (|88p. (|95p . These terms have common origins and the evaluation of their 

limits is essentially the same. The only difference lies in the fact that the shift function, 
F, for the creation/annihilation operators, accounts for an extra particle or hole, as 
described earlier. 

Let us start with the density form factor case 
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We expect the first group of terms to be finite in the TDL because there are no 
obvious singularities in the individual terms. We start by breaking the term into four 
parts as follows 
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(98) 
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where the prime notation is the same as used in the earlier sections, see definitions 
before Eq. (21). Note that because we are only interested in the absolute value when 
evaluating the term Mi, we can freely replace terms by their complex conjugate in 
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order to avoid generating unbounded imaginary parts in the various terms, see e.g. 
p\. The n quasimomenta corresponding to excitations are treated separately and are 
indexed by h. 

We start with the first two pairs of numerator and denominator terms 
• // Xj — /j,k + ic 
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To evaluate the products above we do a logarithmic expansion: 
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(102) 



Similarly, the term p$ can be reduced to the following form 
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Thus, the overall factor M\ takes the following form 
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x exp 
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(X-n h + ic)(X - + ic) 



To obtain the analogous term for boson creation/annihilation form factor, we use the 
same expression as above in Eq. (|104j) . but substitute the shift function F(X), with 
F ± (X). 



4-4- Obtaining term Mi 

Here we will focus on the term 
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This term is expected to contain a power law divergence and a prefactor, and appears in 
all form factors. We require a consistent and physically clear method of extracting the 
divergence and all contributions to the prefactor which are finite in the thermodynamic 
limit. The intuition we will use to do this comes from the fact that we know the 
final answer is expected to scale as a power law of L. In the thermodynamic limit, 
the information about L contained in the quantization of A,/i. Consequently it is 
reasonable to expect that terms which collect to give this power law divergence must 
appear as a difference in Xj — Xk when Xj and Xk are near one another. Therefore 
it is important to isolate such terms on the basis of the nearness of the X's. Thus 
we propose to use a cutoff in order to control the nearness of these parameters and 
perform controlled expansions to extract the divergence and the prefactor. The cutoff 
will drop out of the final answer. This is made more precise below. 

The term appears in the form factors associated with the density operator 
and with the creation/annihilation operators, with a slight modification accounting 
for the difference. We will explicitly calculate the density form factor and point out 
how to modify the result for the creation/annihilation form factors. Furthermore we 
will consider first a single particle-hole pair, with quasimomenta, and indices, 



and lastly generalize the calculation to include multiple excitations. 



It is easier to calculate the thermodynamic limit of A/2 after regrouping terms in 
a convenient way. Our procedure will be the following. 

First we rewrite the terms in round brackets in Eq. (|105|) as T" as shown below 
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The prime notation is the same as used earlier - 
corresponding to particles and holes. 







i.e. we leave off the quasimomenta 
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It is more convenient to evaluate a "temporarily incorrect" term, which we will 
call T", where we include pT instead of p + and then correct the mistake with a 
separate term. Thus to obtain the correct T', we will calculate a correction term T^oie 
that must be multiplied to T" in order to remove the extra terms: 

r = T" x T hole (107) 

The terms that we need to remove are exactly those in which the difference, p^ — A^ 
appears. Consequently, we obtain T^oie in the following way: 
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In the second line, we have relabeled index k as j to obtain line 3 in Eq. (|108l) . 

Finally, the terms that were removed must be replaced by the correct terms. Let 
us call the product of these correct terms T part icie- The way to obtain this term is 
similar to the process used to remove the incorrect terms as in Eq. (|108j) . We get 
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Thus, the final answer will be given by 

M 2 = T" X Thole X T par ticle- 

We start with the expression 
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At this stage there are two levels of approximations that can be made. One, we 
can express differences in A using 

j -k = L {x(Xj) - x(X k )) , 
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Two, the products in Eq. (|llll) above can be exponentiated, and the resulting 
logarithms expanded order by order in 1/L. 

Both these procedures have a region of validity. For instance, it is incorrect to 
perform the logarithmic expansion when j and k are near one another. On the other 
hand, when j and k are far away, the approximation in Eq. (|112[) breaks down. We 
calculate T" by combining these two approximations. The idea behind this is to 
distinguish two regions of the range of Xj, X k using a cutoff - a region when they are 
near one another, and a region when they are far away. Then the two approximations 
outlined above can be used in the region in which they are valid. 

This is illustrated in Fig. OH The blue region I, corresponds to the case when j and 
k are far enough from each other for the logarithmic expansion to become valid. The 
red region II corresponds to when j and k are near enough for Eq. (|112|) to become 
valid. The third, yellow region III will be treated more carefully since Eq. (|112|) is valid 
here, but proximity to the quasi-Fermi points plays a role in how terms are treated 
in this region. Moreover, the denotations a and b are used to differentiate the region 
where j > k and j < k, respectively, to make the calculation more convenient. 

The separation of the regions described above is not arbitrary. To do this we 
introduce a cutoff in the following way - in region I, we only allow the index k to get 
within n* 3> 1 of the index j for all products/sums. In the figure this corresponds to 
the lines at the interface of the red and blue regions. Since we will be working in the 
continuum limit in region I, we would like to know what this condition about indices 
means for the quasi momenta. Since n* is at the crossover where the approximation 
(|112j) just starts to break down, we will use this approximation to define the v* 
corresponding to n* 

j -j± ii* 



v*(X) = \Xj-X k \ 



LplW 



(113) 



Lp L (Xj = A) 

Note that we have a choice of either keeping n* or v* constant - fixing one parameter 
endows the other with a dependence on its position in the range [-q,q] or [1,N] 
respectively. Our particular choice keeps n* fixed, and makes v*{X) dependent on 
A. 
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Using the two approximations in the relevant regions should produce, first, a 
(principal value) integral which captures the contribution to the prefactor from the 
"far-away" region (I a,b in Fig|3]). Second, a constant from the entire region where 
Xj is near Afc, but when they are not both at the same quasi- Fermi point (II a,b in 
Fig. [3]). Lastly, a constant and a power law divergence from the region where both A 
are near the same quasi- Fermi point (III a,b in Fig. [3]). While the intermediate steps 
may produce cutoff dependent terms, these terms should mutually cancel. The final 
answer will be shown to be independent of the cutoff parameter, n* ^> 1. 




Figure 3. Schematic of range of A's. Cutoff parameter n* is used to separate 
regions I. II and III. In region I the indices j, k are far enough to perform 
logarithmic expansions as in, for example, Eq. 11141 . In region II and III, the 
approximation given by Eq. I|112jl is valid, but the logarithmic expansion fails. 
Consequently products must be carefully and explicitly evaluated discretely. In 
region III proximity to the quasi-Fermi point demands greater care in evaluating 
products as demonstrated in Eq. J 1 281) . for example. 



Let us denote the value of T" in region I as Tj. To treat the product in region 
I we may expand the products in T" logarithmically in orders of l/L. Only terms of 
order l/L 2 and higher survive after the double products are evaluated. Without any 
additional approximation, we obtain the following term 

T I= n f i+ . ,f L \ / : k ' . ; r 1/2 u 



F L (X k ) 
Lp L (Xk)(Xj - A fe ) 



F L (X k ) 



L PL {X k ){X 3 - A fc ) L PL {\j){\j - Afc) 



1 - 



Lp L (Xj)(Xj - Afc) 



-1/2 



exp 



1 



Afc) 5 



1/2 



(114) 
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With the definition of the cutoff in A, Eq. (|113[) . we may write 

g{ I] 2^ k L*p L (\ j ) PL (\ k )(\ j -\ k y 




(A-/i) 2 J_ Q Jx+u-{\) 



1 rl r\-V*(\) / i 

-/ Fl(X)d\ d/i h— 

1 J-a+V(-a\ J-q 



q+v*(-q) J-q 

1 rq~""{q) rq / i 

+ 75 / *£(A)dA / 

2 J-g ./A+I/*(A) 

-ir rf A r" ( ^ M f^k^) 

4 ./-g+„*(-g) J-g V A — /X 

-jp*/' ^fftw-^wy, (U5) 

4 7- g A+**(A) V a -m / 

In the last step we have isolated a straightforwardly integrable (lines 3 and 4) and a 
cutoff dependent part (lines 1 and 2). 
Let us first analyze the term 

1 f q /-A-w*(A) II/" 1 / 1 1 

- / F?(A)dA / d/x^ « = 77 / dxF?(qx) , . . 

2 J- q +,'{- q ) J-q (A-M) 2 2 J_ 1+v * { _ q)/q M " >\v*(qx)/q x+1 



l -Fl{-q)\og{qL) - l - 



•v wi^vr «^ 



9+^*(<?) 



2 y_i x + i ' 

(116) 

where we introduced the generalized principal value integrals denoted by P_ , P+ , -P± , 
defined in Appendix B. Note that such integrals can only be well defined when the 
integrand is made dimensionless. This is achieved by mapping the range of integration 
from (-q,q) -> (-1, 1). 

Similarly we may obtain the second cutoff dependent term of Eq. (|115l) . The 
combined answer in region I becomes 

{j rq-v*(q) j rq ~) 

— / d\p L {\)Fl{\) + — dX PL (X)Fl(X)\. (117) 

n* J- q n* J_ q+U * { _ q) J 

In Region II we know that we can no longer expand the log of the products like 
in the previous section because here all terms are O(l) or bigger. We will consider the 
two sets of products in this region, products in Ha, and, lib, separately. In both these 
products we can approximate Xj — X k as in Eq. (|112[) . The second thing is to treat 
terms like Fl(Xj) — Pr,(Afc) with the same precision with which we treat differences in 
A in this region. Thus we need to evaluate the difference in Fl to 0(1/ L) as follows 

Fl(X 3 ) F L {X k ) d [F L (X)\ 
PL{Xj) PL{Xk) dX\p L {X)J 
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Let us define the following notation, Fj — F^Xj), C° = ( jj^jjpj J i where the 



prime denotes that the derivative has been taken with respect to A and evaluated at 
A = Xj. Similarly, define Cj = — ^ 
Thus we have from Eq. (jllll) 




Lp L (\ k ) I I , , Lp L (\ k ) Lp L (\j) I I ^ Lpt(Aj) 



Ajfe Ajfc / y Ajfe 

(^(A^A^V^MA^A^l - Cl/L))^ 
, ; _ ._ n . (Lp £ (A 3 -)A ifc - ^) 1 /2(i PL (A J )A, fe (l - + f»Va 

g g C? - fc)(i + (J - fc ))( 1 + rgCj ~ fc))' 1 !! - C//L) 1 / 2 

3 -=n«+i fc=j-„. (i - fc - i + j. c | (j _ fc) ) 1/2 (i - fc + ( i + i C b (j _fc ))(1 _ C a /i) ) 1/2 (i - cj/l) 1 / 2 

g g (J ^ fc)(l + r^O' - fc ))(! + - k))-Hl - C*/L)^ 

3 =n*+l k=j-n* (j- k - i + ^c| (j _fe)) 1/2 (i - fc + ( i + ^c^ (i _fc)) ( i_c»/L)) 1/2 ( 1 ~ C j/ L ) 1/2 

= n n f — - f 

j=»*+i fc=j-n* U - k - i + icfo'-fe)) 1/2 ^' - fc + (i+icjO'-fc))(i-c;/L)^ 1/2 

g g (J-fc)(l^+^-^)- 1/2 (l~^^)- 1/2 

SL Tjn* + l)r(l - Fj + ±FjC^T(l + F,- + £(gCg + F^))^ 

jJkhi T{n * + 1 - F ^ + r^) 1/2r V + 1 + + W] c ) + ^c?)) 1/2 ' 

In the above expressions, only 0(1/ L) terms will survive and consequently terms 
of higher order in 1/L have been dropped in the intermediate steps. Similarly for the 
region lib 

Tnb = w g g* o--fc)(^+^^)- 1/2 (i-^)- 1/2 

i=i fc=j+i C? - fc - i + iF 3C b) 1/2 (i - fc + i_ i C ai i F c t) 1/2 

= N y( IV + + ^ - {F*C])^T{1 - Fj - \{F 2 C) + FjCfj)^ 

\\ T(n* + 1+Fj- ^F/Cj)Var(n* + 1 - F,- - ^(F^Cj 1 + FjCf))V 2 ' K ' 

We may expand the log of the cutoff dependent terms in the n* ^> 1 limit. This 
leaves the following terms 

Lj^ dAFf (A)pl(A) L p_ q+v , { _ q) dXFj WplW +q (1 



(121) 

The first two terms exactly cancel the cutoff dependent term from region I in 
Eq. (|117l) . The last term can be ignored with the precision to which we are considering 
expansions in the parameter n* . 
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Now we focus on the terms that are not cutoff dependent. First we note that all 
the Cj terms drop out to 0(1/L). We would like to treat the terms. Using the 
notation, ij)(x) = 4-\og(T(x)), we expand the terms to precision 0(1/ L) as follows 



r(i±- 



dx 



)) 



exp[log(r(l± F j ±F j C j /L))} 

exp [log(r(l ± Fj)) ± ± F^FjCj/L] 



= r(l±F L (A J ))exp 



±^{i±f l {x 3 ))\f l {x j ) f^TT^' 

L \PL\Xj) 



(122) 



We will need to evaluate the product over j for such terms. This can be done by first 
noting as a consequence of the chain rule 



^iog(r(x))^io g (r(^A)))x ' 



dx(\) 
dX 



We will evaluate the product of such terms as follows 



(123) 



N-n* 
j=n* 



( ' F > log(rOr))U=i + ^(A 3 .) - £ log(r(x))U =1 _ FL(A3) 



L 



= exp 



exp 



= cxp 



= cxp 



E 



E 



C 'i F 3 



d 



LF' L {\ 3 )d\ 



dx 



{\og{Y{l-F L {Xj))Y{l + F L {Xj))) 



C 3 F 3 







LF' L {\ 3 ) dX 



lop 



Bm{nF L (\j)) 



- q+ J XPL{X) FW)^A^JX))dX 



log 



-(A) \ 



[ 9 d x PLW FLiX)d fFLiX) 

J -q 



d 



f' l (\) d\ \pl(x)J dx 



log 



sin(7rF L (A)) ) 
ttF l (X) \ 



sm(irF L (X))J 



(124) 



Note that we have let v* go to zero. In so doing, we have added in some terms that 
occur at the edges, i.e. form 1 to n* and from N — n* to N. However, these terms 
only make an order 1/L contribution and can be safely treated. Moreover, we treat 
the T(l + Fj)T(l — Fj) terms by extending the product from 1 to N and dividing by 
the "extra terms" 

TV— n* N 

n (m+Fjni-Fj)) 1 / 2 n (m+Fjni-Fj))^ 2 

j=l j=n*+l 



N 

n 



kFl(Xj) \ ( sm(7rF L (q))sm(nF L (-q)) 



sm(irF L (Xj)) 



^F L {~q)F L {q) 



72 



(125) 



where, in the last step we have set all the Fj within n* of the edges to be F^^zLq) since 
this only generates a 1/L error that can be ignored. We have to carry out the above 
procedure since terms involving T(l± Fj) are 0(1) and any terms added or removed 
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may result in a finite constant in the final answer. We combine all the terms obtained 
in region II as follows 



r " - n U(^(A,))j exp b y_/ A ^ (A) n(A) (ma) J ^ v log 



xexp 



— d\p L (\)Fl(\) - - / dA PL (A)F|(A) 

n n J_,j + „*(_<j) 



sin(7rFi(A)) 



]_ \ »*/2 



r(l - F i (g))r(l + F L (q))T(l + F L (-q))T(l - F L {-q)) , 

In this region both Aj, A^ are near the same quasi-Fermi point. Moreover, because 
there are only at most n* terms to consider in this region, we do not require the 
precision we used in the rest of region II. This is because n* / L can be neglected in 
the thermodynamic limit. Thus, by ignoring l/L terms in the near edge region, and, 
simply sending v* to zero in the integral in Eq. (|126[) . we will only be making 0(l/L) 
errors in the final answer. 

This simplifies the analysis considerably, since we will be setting Fj — Fk — 
-Fi(±g) in this region. We obtain two sets of products 

rp _ TT IT 0' - k ) 

HIa f = \ tk ti~ k + F L(-1)) 1/2 (J k Fd-q)) 1 ' 2 ' 

T IIIb = n n — ))1/2( ._ fc _ F n ))1/2 - (127) 

j=N-n*+l k=j+l w L-Vi// \J l-VI)) 

Starting with the first expression and expanding the answer for n* 3> 1, we obtain 

rp = rq)r(i + F L {-g))^T{i - FU-g)) 1 / 2 

IIIa f = \ nj + F L (-q)y/ 2 r(j - F L (- q )y/ 2 

= G(n* + 1)G(1 + F L (-g)) 1 / 2 G(I - F L (~g)) 1 / 2 r(I + F L (-g)) n ' / 2 T(1 - F L (-g))»V 2 
G{n* + F L (-q) + l)V 2 G(n* - F L (-q) + 1)V 2 

« T(l + F L (-q)) n ' /2 r(l - F L {-q)) n '/ 2 G{l - F L (-q))^ 2 G(l + ^(-q)) 1 / 2 ^*)-^, 

(128) 

where G{x) denotes the Barnes G function, (see Eqs.(Al) - (A4)). 
Similarly, for the other product 

T im = r(l + F L {q)Y'l 2 T{l - F L {q)) n */ 2 G(l - F L (q)) 1/2 G(l + F L (q))^ 2 (n*)-^. 



(126) 



(129) 

The final result for T" is given by 
T" = (qLp L (q))-^ F ^ )+F - { - q)) (G(l + F L (q))G(l + F L (-q))G(l - F L (q))G(l - F^-q))) 1 / 2 



X 



JJl^FL^i))) \ J -i x 2 ~l Aj_ q J_ q 1 \ A- M J 
\ fd\{l ~ ^(A)cot (7 r^(A))} [f' l {X) - Fl[ ^ X) ) 



(130) 
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Having obtained the expression for T", we need to remove the "incorrect" hole 
terms, Eq. (I108|) . to obtain T' as in Eq. (| 107|) . The treatment of this term depends 
on the distance of the hole from the quasi-Fermi points. To quantify the notion of the 
distance from the nearest quasi-Fermi point, we define a quantity qk in the following 

way 

q~l = % - n N - 1. (131) 

Thus if the k th hole is at the right quasi-Fermi point, qk = — 1, if it is in the middle 
of distribution then qk ~ —N/2 etc. Moreover, we say that a hole is near the right 
quasi-Fermi point, if qu/L is in the thermodynamic limit, and conversely we say 
it is "deep", or far from the right quasi-Fermi point, if q^/L is a constant in the 
thermodynamic limit. 

Similarly for holes near the left quasi-Fermi point, we define q^ to be 

% = n k ~ n i + 1- (132) 
Now we wish to evaluate 



T ho i e = n 1 



(133) 



We would again like to use the cutoff procedure used to obtain T". In the region 
where \ j — i~\ > n* , we expand Thole logarithmically as follows 



T hole = CX P 



E lQ g 



1 



Fl{X 3 



\j-i-\>n* 



L PL (X 3 )(X 3 -X-)+ F ^-J^ 



exp 



E 



Fl{X 3 ) 



r L PL {\ 3 ){\ 3 -\-) 



O 



L 



(134) 



Jj-i \>n* 

Meanwhile, in the region where \j — i~\ < n*, we use the approximation in 
Eq. (|112[) . to evalute discrete products as follows 

F(\- 



rpll _ 

hole 



n 



J 



(135) 



Case I : Hole is deep inside the distribution 

When the hole is deep inside the distribution, we obtain from Eqs. (I134[) . (jl35[) 



log(T, 



k,deep\ 

hole / 



CI A 

x-x- 



-9 



oixl^ 



log 



T(l + F L (\7))T(l-F L (X-))T(n*) 2 ) 



(136) 



T(n* + F L (X^))T(n* - F L (X^)) 

We get cancellation between divergent and cutoff dependent terms from the integrals 
due to the symmetric cutoff. Thus in the thermodynamic limit we have 



,deep _ sin[7rF(/^ fc )] 



" hole 



exp P 



dx 



F{qx) 



Case II: Hole near the quasi-Fermi point 



(137) 
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Meanwhile when the hole is near either ±g, we obtain from Eqs. (I134p . (|135p . the 
following expression 



Tk , ed9 e (aLo(a)) TF(± q ) MnF(±q)}r( T q±±F(±q)) f f 1 F(qx) 

where we use to refer to the location of holes near a quasi-Fermi point as defined 
in Eqs. JEEJ, ([T52]> . 

In order to replace the "correct" term in the evaluation of M2, we require the 
following term: 



Tparticie 



(139) 



Again, we distinguish the case where the particle excitation is far from either edge and 
the case where it is near the edge. For this purpose we define a means of quantifying 
nearness to a quasi-Fermi point. Let us define, = — njv — 1, for a particle 
near the right quasi-Fermi point, and p^ = — ni + 1, for a particle near the left 
quasi-Fermi point. We treat two distinct cases, 
Case I: Particle far from edge 

If we find Pk/L is finite in the thermodynamic limit, then we have the following 
convergent term 

(140) 




T 



particle 



Case II: Particle near quasi-Fermi Point 

On the other hand if P^/L vanishes in the thermodynamic limit, we need to 
use our cutoff procedure to separately evaluate the contribution to T partic i e from the 
region where A^ is near p + , and where they are far. We let the j product extend as 
near as v* (±g) to the relevant quasi-Fermi point, and treat the terms within v* of the 
quasi-Fermi point separately 

, = / qLp(q) \ ±Fi±q) T(n*)r(±p± + 1 T F(±q)) 
partlcle \ n* J r(±p± + l)r(n* TF (±q)) 

= (qL P (q)f^ r ^ + 1 /^»ex P f-P ± f 1 dx^M) . 

1 Iy " r(± P ± + 1) y \ ± J_ 1 xTi) 

Now we consider how the term M2 needs to be modified when there are n particle- 
hole pairs. Most obviously, the shift function used, ^(A), becomes the sum of the shift 
functions with the different particle-hole contributions. The term T" is unaltered since 
it is evaluated with substituted for pf , making the definition of F(X) in Eq. ([7^)1 
valid for all Xj. In order to "correct" for this convenient substitution, we include 
the terms T^oie associated with each of the holes, and terms T part icie associated with 
each of the particles as described in the previous subsections. However, when there 
are multiple particle- hole pairs, this is no longer sufficient to correctly evaluate Mi- 
There will be another group of terms which we shall call T cross whose origin can be 
understood in the following way. Let us start by writing down an expression for Mi 
and show how to arrive at the correct final answer 



(141) 



M 2 = T"xl[T^ Ucle xTi k l 

k=\ 
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( N 




U 




U 





= T " x ( II TpaLde x nlj x T cross , (142) 

where the notation T^] e (T^l ticle ) means the contribution from the hole (particle) 
term specific to the details of the excitation, and refers to the index 

corresponding to the excitation pair /i^ , ■ 
Thus we obtain the expression 

Cr0SS = ll\ K-% + J 2 ) ■ (43) 

We expect that the contributions to T part i C i e , T^oie from the excitations 
{/i^ , fall into one of the cases discussed in previous subsections. 

Thus, the form of these terms is always known in principle. 

The term T cross is sensitive to the details of the excitations - for instance it will 
contain divergences when some number of the particles, or holes, are clustered near 
each other. For our purposes, this is only relevant to form factors of states containing 
high order Umklapp excitations, i.e. several particles and holes near the left (right) 
quasi-Fermi points. We will discuss the term T cross in more detail in the section 
devoted to Umklapp form factors. 

The discussion above can be made applicable to Mi in the creation/annihilation 
operators with two modifications. First we need to replace Fl(X) by P± l(A) defined 
in Eqs. (|89)l . (|90|) . as discussed in the section following these definitions. 

Secondly, if we wish to retain the calculation performed in the previous subsection 
then we need to extend/reduce the upper bounds of the indices of the products in M2 
(see Eq. (fTOSjl ) from N-l to N by using the relation, fi~ +1 = X n - - l p ^}) +0(1/ L 2 ) 
to define the "missing" ground state or excited state quasi momentum (see Section 
IV. 3), and for concreteness we again concentrate on the case of the annihilation 
operator. This will introduce some error to be corrected to obtain the final answer for 
M2 and - for the creation and annihilation operator respectively. 

We will follow the consequences of such a replacement and correction for the 
annihilation operator and just present the final answer for the creation operator since 
the derivation is essentially the same. 

Consider the product 

1/2 
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Now extending the upper limit of indices to N introduces the following extra terms 



N 

n 



(A n - - Aj)(/x, 



n+l l l j) 



(A n - -frj){X 3 -m„+i) 



(145) 



Notice, however that one group of terms in the above expression is present in the final 
answer for the annihilation operator, , which can be seen explicitly in Eq. (|93[) . 
Thus we may write 

1/2 
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An- — Aj 



Xr, 
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AT 
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(Xj - fi k ) 2 
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n 



Mn+l 



-A, 



(146) 



where the last product serves to cancel the unwanted terms in the parentheses. For 
notational convenience let us call the product multiplying on the left hand side of 
the above equation, Tf xtra . However, as before it is more convenient to evaluate our 
double primed products where we substitute /i^ in place of /if and correct this mistake 
separately as in Eq. (|107p . Consequently we obtain from Eq. (|146p and Eq. (|142[) 
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(147) 
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fci^fc 2 



n+l n 
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(148) 



n 

fci^fc 2 fc 2 =i 

The corresponding term for the creation operator form factor can be obtained by 
switching particle and hole quasimomenta in the above expression. 

We obtain a similar expression for the creation operator and the answers are 
summarized in the next section. 

The most general result for M% in the case of the density form factor is given by 

M 2 = {qLp L {q))-^ F ^ +F ^-^ (G(l + F L (q))G(l + F L (-q))G(l - F L (q))G(l - F L {-q))f' 2 



nF^Xj 



xexp 
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hole particle 



(149) 



Similarly, for the creation/annihilation operators it is given by 



(qLp L (q)) 



(G(l + F ±tL {q))G{\ + F ± , L (-q))G(l - F ±>L (q))G(l - F ±>L (-q))) 
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xexp 



\ f q d\{l - nF ±>L (X)cot(7rF ± , L (m (^W - ^j^f^' 

n(+l) ra+l(-l) 

n T W y TT T ( fe ) 
J-hole A 11 1 particle 

k=l k=l 



Tcrots ■ (150) 



^.5. Fredholm Determinants 

The last terms left to express in the thermodynamic limit are the determinant terms 
appearing in the form factors, such as 9 in Eq. ([86)) . We would like to write 
such terms as Fredholm determinants (see Appendix D and Ref. [75 for details on 
Fredholm determinants). The determinants are slightly different for the density and 
creation/annihilation operators. We will first do the calculation for the determinant in 
the density form factor, and then use some of these results to express the determinant 
for the creation/annihilation form factors. 
We need to evaluate 

Det (s jk + EL^ ^(K(Xj, A fc ) - K(X P , X k )) 

e = — ^ 3 - — 3 - T . 

v+ - v p - 

(151) 

Due to properties of the matrices involved in the determinant, the choice of X p is 
entirely arbitrary [46] and it need not be from the set {Xi}. We can use this fact 
by taking X p to be much larger than all other parameters of the problem. Under 

these circumstances, K(X P , Xk) — > f§- Meanwhile the denominator lends itself to the 

p 

following expansion 

y± = TT (l - M * ~ A * ( X h -vt)( X h ±ic Y 



TT" h a F ( A ") i F{X m ){X m ±ic) \ (I \ 
n ( l+ LAX^X p + L P (X m )Xp J + °UJ (ld2) 



Using this, one can evaluate the denominator in this limit as 



V p + -V p = 2 ^ [Y,K-ti+ I* F(X)dx)=^f^, (153) 

where we have used the relation between momentum and quasimomenta as in Ref. [43. . 

The numerator also lends itself to a convenient re-expression. Essentially the 
numerator has the form Det(A + a—), where 

_ i(nj - Xj) -pi- \i m - Aj 

2c 2 

a =9-. (154) 

Using the fact that B is a matrix of rank one, we can expand the determinant using 
rows and minors alternately from A and B. Using the expression for the determinant 
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of a sum of matrices [43] . we find that because of -B's rank, only terms with one row 
from B will survive. Moreover the determinant of A also vanishes permitting further 
simplification. Thus the determinant can be written 

Det(A + aB) = Det(A) + a B^a^ + a Ay ay — a Ay-ay 
= Det(A) - aDet(A) + aDet(A + B) 

= aDet(A + B). (155) 
Taken together with the previous statements we get 

e - (* + fff n b^wr* - 1)) ■ ™ 

To express the above as a Fredholm determinant, we consider the following term 

The terms V^- admit the following partial limit 

J M - A, ± zc eXP A - Ai ± 

The infinite product in Eq. (|156[) is not as easy to treat. We cannot make the same 
kind of continuum approximation as above because of the bad behavior of such an 
integral around Xj . This can be treated by temporarily introducing a cutoff into the 
continuum version of the product (an integral) and extracting any constant (cutoff 
independent) error committed in so doing to finally obtain an answer that is cutoff 
independent. We obtain 

TT ( Mm - Aj \ _ TT 1 ( \ — F UKn) \ f /i + - A 3 



0(1/L) 



(158) 



exp 



' _ f q l r(n* + l + Fj)iy + l-F 3 -) 

, a - a, j Xr x-xA r(i + F,-)r(n* + i)i>* + i)r(i - Fj) 



X 



n / + \ 
TT [ ^fe ~ A 3 



(159) 



Here Al and A# correspond to the spectral parameters corresponding to the quantum 
number j ± n*. Under the assumption that n* 3> 1, we can use the Stirling 
approximation on the gamma functions: 

r(n* + l + Fj)r(n* + l-Fj) 

r(i + F i )r(n* + i)r(n* + i)r(i-F i ) 

= exp {(n* + 1 + Fj)log(n* + 1 + Fj) - n - 1 - F 5 + (n* + 1 - Fj)log{n* + 1 - F,)} 
x exp {-„-! + Fj - 2(n* + + D + 2(» * +D> r(l + F £ (Aj))r(l - F L (Xj)) 



1 
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Where the exponential terms exactly cancel out. Furthermore the integrals contribute 
terms at the boundaries that also mutually cancel, leaving no constant terms upto 
0(1/L) 

houndaries = - F(Xj)log(X L - Xj) + F(Xj)log(X R - Xj) 



a, = 



= 0. 

This allows us to express aj as 
i(fij - Xj) 



log 



Lp{X 3 ) 



log 



Lp(Xj) 



(161) 



nLr 



A 7 —\j+ic 



exp 



r dx- 



n 



/4 - Aj 



exp 



Pf 1 dx zm 

j — i x ~^j 



/'i 



n \- X 

IT 



-g ' A— Aj -\-ic 

ic Mft - A 



C.C h = 



r(i + F(Aj))r(i — F(Xj)) 



X 



2nLp{Xj) M + - Aj + ic A h - Aj 



x exp 



-P 



F( ga Q , 

CtX ; h Z7T-T (A,* 

-1 x-Xj/q 



F(X) 



dX^ x 

A A, • /<• 



Let us define an operator, G that acts on [— q, q] x [—5, g] and is given by 
GO, v) = a(/x) \K(fi - v) - - J , 



where the function a(/i) appearing above is 



a(M) = 



l n \- 

TT — 

h=i Mfc - H + ic X h -n 



2tt 



ri + ic uT — a 

x — exp 



-P 



F{qx) 



dx- 
1 x- fi/q 



iTrF(^) 



(162) 
(163) 

V ^ 

q X - [i + ic 



(164) 



Using Eqs. ([T56]l . (fl63]l . we obtain 
6 = -^Det (/ + g). 

(165) 

The determinant terms appearing in the creation/annihilation form factors are a 
slightly modified version of the one for the density form factor, and are given by 



(A? ~ Aj) Hm^j^n- Aj) 



C" 



Det S 



V7 - V7 



rim^i CAni - Aj 



(X(Aj - A fc ) - K{n 



n+1 



A fe )) , 



(+ = Det *, fc 



jjfij - Aj) nm^jtAro - Aj) 
Vj -V j llm^j(Am-Aj) 



(K(X 3 - A*) - if (Aj - M + +1 )) 



(166) 



where the +/— in the subscript correspond to the determinant for the 

= ± 

creation/annihilation operators respectively. Note that by Vj we mean the term 
occurring in the creation operator form factor defined analogously to Eq. (|57j) . 

The procedure we use to express as a Fredholm determinant is similar to what 
we used for the case of the density form factor. Let us define functions similar to 
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i(fij - Xj) rim^iCAn 



A,) 



N 

N-lf~ 



vf-vrnZ^Xm-Xjy 



- Xj) Tlm-Lj {Prn - Aj) 



= + 



V n m /j(Ar, 



A,) 



(167) 



In the thermodynamic limit, and to the same accuracy with which we calculated 
a(/i) we find 



n+1 . _ . n 



x exp 



-P 



F-(qx) 



dx- 
1 x- fi/q 



mF_ (/i) 



9 dA- 



X — fi + ic 



6+( M ) 



1 A A7 



n 



fx + ic 



x exp 



-P 



i x-fi/q 



Ii 



H + ic 



iirF + {n) 



q dX- 



A — n + ic 

(168) 

Now we may evaluate as Fredholm determinants of I + H ± {^ 1 v), where 
act on (—q,q) and are given by 

H-(jm, v) = b-((j,)(K(n -v)- K(»- +1 - i/)), 

fr+(/i, !/) - b + {n){K{n - u) - K{ll - n+)). (169) 



4-6. Thermodynamic Limit of Form Factors 

We now present form factors for a few excited states that we would like to relate to the 
correlation prefactors of the Lieb-Liniger model. We first present a general expression 
for the density form factor and then obtain specific form factors from it, and then 
repeat the procedure with the creation/annihilation form factors. 

By combining the results from Eqs. ([88]). {EH}, ([139]) . ([IBS]) , and the 
simplification, (|C.3[) . described in Appendix C, we obtain a general expression for 
the density form factor: 

T = c {qLp{q))-^ F2{q)+F2{ - q)) (27r)^)- F (-«»G(l - F(q))G(l + F(-q)) 



exp < P ± 



dx 



F 2 (qx) 



exp 



exp 



E 

i=l 
1 

.fe=i 



-i 

g 
-q 

8 



.x 2 - 1 



r/A 



-9 

f(a)( m +-mD 



V f(a) - f(m) 

-Q 



A — /i 



(A +ic)(A-/i^+ic) 



n 



-^ + +»c) 2 

OiTj + ic )it4,j + ic ) 



1/2 



d/i / dA 



(*0 T (fc) 

particle hole 



F(X)F( Ll ) 
(A — /i + «c) : 

n 



Det(l + G) 



Dct 



L = -L Vi(p(^)p(^)) 1/2 (M t " 



(170) 
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In the above expressions, L is the length of the system, q the quasi momentum at the 
edge of the distribution, are the quasi momenta corresponding to excitations, 

and, p(X) is the ground state distribution function defined in Eq. ([50]) . Moreover, 
F(X) is the composite shift function defined in Eqs. (I5T|) . (|75|) . G(x) is the Barnes 
function (Appendix A), while G is obtained from Sec. 4.5, K(X) = c -i 2 + x i and the 
terms T par ti C i e , Thole, T cross are described in Sec. 4.4 and cannot be written down in 
a general form. 

Let us first calculate the density form factor for the Umklapp state. This state 
contains m adjacent holes starting at the left quasi-Fermi point, and m particles 
starting at the first available spot after the right quasi-Fermi point. Consequently we 
obtain m contributions from T part icie, as in Eq. (|14ip . of the following form: 



(171) 



n" 1 ^(k) , T , ,,mF(q) jr - F(q)) ( f 1 mF(qx)\ 
T particle = [[ exp -P + / dx I , 

k=l fe=l V ; V J 1 7 

and similar contribution from the product of the T^oie terms, given by Eq. (|138[) . 
TT T (fc) - (aLo(a)) mF{ - q) TT T(k - F(-g)) ( f 1 mF(gx) \ 

(172) 

Next from the cross terms given by Eq. (|143[) . T crossl and using the approximation 
Eq. (|112[) combined with the fact that the holes and particles occupy m adjacent spots 
near the two quasi-Fermi points, we obtain up to 0(1/ L) terms of the following form: 

m m 

n(^) i/2 -(^))- i/2(m2 - m) n r W' 

j^k i=l 
m m 

no4) i/2 =(^))- i/2(m2 - m) n r ( j )' 
&k <=i 

Y[(p+-p-) = (2qr 2 - m , 

T cross = (2Lqp(q))- m2+m G 2 (m + 1). (173) 
Lastly, there is a term coming from the product (for instance appearing in 

ft ( *m ^ = (^my m (174) 

where, we have only retained terms bigger than 0(1/ L) by assuming all the holes and 
particles occur at approximately the quasi-Fermi points. Furthermore we have used 
the symmetry of p(X). 

Combining Eqs. (|171til73|) and using the symmetry of the shift function for the 
Umklapp, and, the relation (Al) we obtain 

rp ( T i \\2mF( 9 )-m 2 G(m + 1 - F(q)) 2 

J Umklapp - (q^P(Q)) ( 2 )m 2 G(l - F ( q ))*T (F (q))™T (1 - F(q))™ 

x exp (p_ f dx^M P+ J" feggM) . (175) 
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Furthermore we may also relate the shift function F(X) to the ground state density 
function p(X) for the Umklapp [5B]. Let us start from the equation for the shift due 
to a first order Umklapp, F{X\ — q) — F(X\q), i.e. the shift for a particle hole pair at 

F(\)--]_toK{\-v)F{y) = — — , 

F(X) -7^f d»K{\ - v)F{v) = ~ J" dvK{\ - u), 

l r q 

F W - 7T / dvK{\ - v){F{v) - 1) = 0. 

27r J-q 

(176) 

But this last equation is solved by F(X) = — (27rp(A) — 1), since 

1 f q 

2tt P (X) - — / dvK(X - v)2tt P {v) -1=0. (177) 

27r J-q 

Moreover, since an order m Umklapp can be constructed by repeatedly performing m 
order one Umklapps, we may use the linearity of F(X) to conclude 

F m (A|Umklapp) = m(F(X\ - q) - F(X\q)) = -m{2np(X) - 1). (178) 

This makes the exponent of L 

— F 2 (q) + 2mF(q) — m 2 = — (m\f~K — m) 2 — 2m(m\f~K — m) — m 2 

= -m 2 K, (179) 

where we have used the relation 

2irp{q) = VK, (180) 

from Ref. [43], where K is the Luttinger parameter [5]. 

We may now obtain A m the coefficients in the equal time density-density 
correlator Eq. ((T|) using Eq. (fTo) . We substitute the various expressions obtained 
above into Eq. (11701) to obtain 

'4<7 2 + c 2 Y"V (7(1 + mVK) 2 G{m + 1 - my/K ) 



A, 




Ac 2 J \ r(m - m\f~K) m Y{l - m + m\f~K) m G(\ - m + m-fK) 
eXp4m9/ d \(X + icr-q*) -J_ q d »J_ q d \x-» + icy 



F 2 {qx) - 2mF{qx) 1 f q f q „ , / F(X) - F(p) 

dXdp 



1 2 J J \ X - ii 



Bet 2 (l + G) 



Det 2 <i K 



2tt 



To determine the prefactors of the singularities of the density structure factor 
(DSF) S{k,uj) at Lieb's colective modes ei, 2 (fe) [38l|43] using Eq. J27]) we need to find 
corresponding form factors of the density operator from Eq. (|170|) . 

By considering the form factor for a Bose gas with one hole at the right edge and 
a high momentum particle, we may determine the prefactors of the DSF singularity 
at £\{k). To do this we need to relate the quasi momentum of the particle /i + to the 
momentum k of the excited bosonic state as follows [43| [21] : 

k = p + - npa - { dX6(p + - X)p(X), (182) 

J-q 



Exact prefactors of correlation functions of ID quantum integrable models 



47 



where, 0(A) = i In I ) , and q is the quasimomentum at the edge of the distribution. 



For the excited state with the particle hole pair (/i + , q), we substitute F(X\p, + , q) 
from Eq. (juTj) and substitute this function in the relevant places in Eq. (1170[) . For 
such a particle-hole pair, we obtain a form for T^oie as in Eq. (|138[> . and for T part i C ie 
as in Eq. (|139p . Since there is only a single particle-hole pair there is no contribution 



S\{k) 




-(F 2 (q)+F 2 (-q)-2F(q)+2) 



exp i P± 



exp 



exp 



2F 2 {qx) 
x 2 - 1 " 

2F(X)(p+-q) 



dx 
-l x 



^F(q)-F(-q) 

2F{qx) 



G(l + F(-q))G(l-F(q)) 



(A- q + ic)(X- p+ + ic) 

q dX- F ^ 



1 

+ ^2 | „2 



dX 



dfi 



T(l-F(q)) 
F(X)-F(jmY 



X — jJL 



(q-H+Y + c 



(A - fi + icf 



dX 



2F(X) 



F 2 (q) 



p(p+)(p+/q-l) 2 



X 



+ 



Det 2 (l + G) 



Dct" 



l- £- 



(183) 



Note that the momentum dependence in the above expression is contained in the 
shift function, F(X) — F(X\p + ) — F(X\q), as well as the terms directly involving p + . 
Furthermore we have used the relation in Eq. (|180[) to express the answer in terms of 
the Luttinger parameter, K [S]. 

Similarly, we may use Eq. (I170p . to obtain the DSF at the Lieb mode £2(^)7 by 
considering the density form factor of a system with a single high momentum hole 
and a low momentum particle (at the right quasi- Fermi point). This time we relate 
the quasi momentum of the hole, \i~ , to the momentum, k using 

(184) 



— k = fi — npo — / dX9(p, — X)p(X), 

J-q 

and, use the unique solution of the above to determine a shift function, F(X\q, p~), 
using Eq. (|5T1) . We may then use this F(X) in Eq. (|170[) in conjunction with the 
expression for T part i C ie given by Eq. (jl41[) . and Thole given by Eq. (jl 36[) to get the 
following expression for S*2(fc): 



S 2 (k) 



-(F 2 (q)+F 2 (-q)+2F(q)+2) 



(o ^(,)-F(-„) ( G{1 + F{-q))G{\ - F{q))T{\ - F(q)) Y 

[ ! { r(i + F(p-))r(i-F(p-)) J 



exp 



exp 



exp 



dX- 



2F(X)(q- f x- 



(A — /x + ic) (A - q + ic) 



(q-^-) 2 + c 2 



F 2 (p~ 



p(^-)(p-/q-l) 2 



dp. 



q dx f (« 



P 



1 dx 2F ^ X ] 

-l x-fi fq 



(A - p + ic) 2 



dX 



dp 



dxK M -? + 



( F(X)-F(p) 
V X - p 



1 dsH^Ml 
1 x-1 J 

Det 2 (l + G) 



Det" 



('-*) 



(185) 



Here too the dependence on momentum, k, is carried by F(X) — F(X\q) — F(X\p~) 
and p~ . Furthermore we have used the relation in Eq. (I180p to express the answer in 
terms of the Luttinger parameter, K [S]. 

Using the results from Eqs. (|M]l. (pp]) . (|104l) . f)150[) as well as the simplification 
(|C.3[) outlined in Appendix C, we obtain the following analytic expressions for the 
form factors of the creation and annihilation operators. 
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The most general expression for the annihilation operator form factor is given by 



F_ 



n 



F-{lh) 



n+l 

nxT 
} 

k=l k=l 



n T (k) 
P< 



(k) 
hole 



x T C1 



exp 



dx 



Fl (qx) 
x 2 - 1 



dx /V F - (A) - F - (/i) 

-g J-q 



A — /i 



exp 



exp 



E 



dA- 



(A — + ic) (A - + ic) 



n 



-,1/2 



F_(A)F_( M ) 
2(A - p + icf 



Det(l + ff-) 



ic )(l4,j 



Dct 1 - 



2tt y 



The creation operator form factor is given by 



g + = 



(<lLp(q)) 



L 1 /2 p ( Ai +)l/2 



n+l 



(fc) 

particle 



exp <( P± / dx^±— T 



.fc=i 



fc=i 



x T„ 



F + (A)-F+0i) 



exp 



exp 



E 



dA- 



F + (A)(/i+-/V) 



,i=l "-9 
■'I 



(A-/i, t + ic)(A-^++ic) 



n 



-9 



2(A- y u + ic) 2 



J J J L v »iJ 

Det(l + H 



A — /i 

{ih,, +ic)(t4j + ic ) 



"I 1/2 



Det 1 



(187) 



In the above expressions, L is the length of the system, q the quasi momentum at the 
edge of the distribution, p~ , p + are the quasi momenta corresponding to excitations, 
and p(X) is the ground state distribution function defined in Eq. (|60l) . Moreover, 
F±(X) are the modified shift functions defined in Eqs. ([TO]) , G(x) is the Barnes 

and the terms 



2 c 



function (Appendix A), is obtained from Sec. 4.6, K(\) = 

Tparticie, Thole, T cross are described in Sec. 4.4 and cannot be written down in a general 
form. 

In order to obtain from Eq. (|186j) the annihilation operator form factor for an 
order m Umklapp state, we first collect terms from Eqs. (|171[) - (|173p with the 
appropriate modified shift function and the extra hole to obtain 



Umklapp 

= (Lp(q)) 



(m-l)F-(q)+mF-(-q)-m(m-l) 



1 



(2?)' 



G(m + 1 - F-(-q))G(m + 1 - F_(g)) 



G(l - F-{q))G(l - F_(- g ))r(F_(-g))™- 1 r(l - F.(-q))^T{l - F_(g))r(l + F_(g)) 



x exp 



-(m- 1)P+ / 

'-l ■'' 



F_( g x) 



f 1 , F_(qx) 
mP- j dx- ' 



x ■ 



1 



(188) 
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For this excitation we may determine the shift function at the quasi-Fermi points, 
F_{±q), using Eq. flU: 

F_{±q) = F,4±g|Umklapp) - F(±q\q). (189) 

We can substitute F m (±<7|Umklapp) = -m(\ r K - 1), from Eq. ([178]) and Eq. (fl~80]) . 
and substitue F (±q\q) + tt p{±q) = 1/2 ± (1/2 - l/(2\/K) from Ref. 21, to determine 
the exponent of L in the Umklapp annihilation form factor, 

- ■xiFliq) + Fl(-q) - 2(m-l)F-(q) - 2mF-{-q) + 2m 2 - 2m + 1) = -m 2 K - —. 
2 AK 

(190) 

where, K is the Luttinger parameter [5]. 

We may then use Eq. (fl~4"|) to relate the absolute square of the annihilation 
operator form factor Q~ for the Umklapp state to determine the coefficients B m of 
the Green's function of the Bose gas in Eq. © to obtain 

y K) - m >K- y *K / V + c^ \ (m ~ 1)2 Det 2 (l + g-) 



X 



exp 



exp 



2^-2m-i v , \ c 2 J Det 2 (l-§ 

G{m - F-(g))G{m + 1 - F_(-q)))G(l + F-(-g)) N 2 
r(F_(-g))™-ir(l - 1G(1 - 

y_ 9 ((A + ic) 2 -q 2 ) J_ q J_ q (X-p + ic) 1 

1 F 2 M-2(m-l/2)^_( ga ;)+xF_(gx) 1 f /"« / f_(A) - F-(p) ^ J 

2P ±/ ^ -JJ_dXd^ X — fi 



(191) 



Note that we have used the relation in Eq. (|180[) to express the answer in terms 
of the Luttinger parameter, K [5]. 

To determine the prefactors of the singularities of the spectral function A(k,uS) 
near Lieb's colective modes ±ei.2(fc) [MISS] we need to find corresponding form factors 
of the creation/annihilation operators from Eqs. (|186[) . (|187[) . 

We obtain A+(k) from the creation operator form factor for a Bose gas with a 
high energy particle, and obtain A- (k) from the annihilation operator form factor 
for a Bose gas with one high energy hole. These functions give the prefactors for the 
spectral function singularities at ei(fc) and — 62(h) respectively. We must first relate 
the quasi momentum of the particle to the momentum k of the excited bosonic state 
as follows g3[H]: 

k = p + - upa - f 9 dXd(p + - X)p(X), (192) 

and for the quasi momentum of the hole 

- k = p- - irpo - ( dX9(p- - A)p(A), (193) 

J-q 



where 9(X) = i ln^ ^-a ) ' an< ^ 1 ^ s ^ ne Q[uasimomenturja at the edge of the distribution. 

Particle-hole pairs in Eqs. (|8"9")l . (|9"0"|) are defined with respect to ground states 
of N ± 1 particles, so to obtain a state with a single high energy particle (hole), a 
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particle-hole pair needs to contain a hole (particle) at q. Such a procedure results in 
the following expressions for shift functions F± (A) 

F^(Xj = F(X\p + ) + Trp(X), 

F_(X) = —F(X\fi~) - np(X). (194) 
We may then substitute these functions in the relevant places in Eq. (|187p and 
Eq. (I186p . combine terms from Tj lo i e given by Eq. (| 138[) . and T part i C ie given by 
Eq. fI59"]) . to obtain 



A±(k) = 



q jK)- (Fl[q)+Fl{ - q)) { 2,f^- F ^-\G{l + F^))G(l-FW))f 
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exp 



exp 



dp \ dX 
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dX- 



(X - p + ic) 2 



Det z (l +H-) {{q-^f + c 2 ) 



Det' 



(A — /i* + ic)(X - q + ic) 

Note that the momentum dependence in the above expressions is contained in the shift 
function, -F±(A) defined above, as well as the terms directly involving \x . Furthermore 
we have used the relation in Eq. (|180p to express the answer in terms of the Luttinger 
parameter, K [5]. 

We obtain the prefactors of the spectral function singularities, A+(k) and A-(k), 
near the Licb modes £2(fc) and — ei(fc) from the form factors of the creation and 
annihilation operator respectively. We can obtain the pref actor A+lk) from the 
creation operator form factor of a system with a high momentum hole and two particles 
at the right quasi-Fermi point. From the annihilation operator form factor of a system 
with a high momentum particle and two holes at the right quasi-Fermi point, we obtain 
the prefactor A_ (k). We again relate the quasimomenta of the particle and hole to the 
momentum k using, Eq. (|192j) and Eq. (|193[) and determine F±(X) using Eqs. ([89]l. 

Particle-hole pairs in Eqs. (l89|) . (|90|) arc defined with respect to ground states 
of N ± 1 particles, so to obtain a state with a single high energy particle (hole), a 
particle-hole pair needs to contain a hole (particle) at q. Such a procedure results in 
the following expressions for shift functions F±(X): 

2F(X\q) 

F(X\^ 



F + (X) 



F_(A) 



") 



F(A| M -)+7rp(A), 
-2F(A|g)-7rp(A). 



(196) 



Substituting these results in Eq. (|187j) and Eq. (|186|) in conjunction with the 
expression for T part i C ie given by Eq. (I14ip . and Thole given by Eq. f| 136[) we get the 
following expressions: 



A+(k) 



p{n-)(v-/q-iy 

fG(l+F + (-q))G(l 



-(FUq)+FU-q)-4F + (q)+4) 



(195) 



(27* 



-(g)-F + (-g)-2((g-M ) 2 + C 2 ) 



F + (q))T(l-F + (q))T(2-F + (q)y 



Det 2 {l + H+) 



r(i - F + (/i-))r(F + (M-)) 



Det 



(!-*) 
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F + (X)(q-^-) 



"i 2Fl(qx) 
exp { P± I dx = — — P + 



-l 



1 



1 4 F + ( g x) /•« 
cfe ; - + 2 dX 



-l 



a: - 1 



(A — fi + ic) (A — q + ic) 



exp 



exp 



1 2F+(qx) 
P I dX = 

l x-fj, /q 



-q J-q y 



A — /i 



9 d r dA F+(A) F+(») 

-a J -a (A - M + «C) 2 



(197) 



-(F* (g)+F=(-«)+4F_(g)+4) 



(2tt 



(-S)-2 



p( / i+)( A i+/ ( ?-l)4 

Dct 2 (l + H~) ((q - /i+) 2 + c 2 ) / G(l + JL(-g) )G(l - F_(g) )r(l + F-(g) )r(2 + K- 



(«))' 



Dot' 



r(F_( 9 ))r(l-F_(g)) 



r-i 2Fl{qx) 



1 4F_( ? a:) 
dx- 

-l 



x- 1 



-9 ^-9 V 



X — fl 



exp 



exp 



9 , r q F-M F -(M) 

dfi / aA - 

-g ./-g (A-^i + Zc) 2 

2/ ■-, F -(^ + -^ 

q (X — q + ic)(X — \i + + ic) 



i 2F_ (A) 
dX ^T 

-a X - fl+ 



Here too the dependence on momentum, k, is carried by F±(X) defined above 
and fjr 1 . Furthermore we have used the relation in Eq. (|180[) to express the answer in 
terms of the Luttinger parameter, K [S]. 



(198) 



4-7. Numerical Results 

Let us plot here some prefactors which are obtained using above analytical results. 

A few prefactors of equal time correlators A m , B m , see Eqs. ([1} and ([2]) are plotted 
in Fig. U] 

Si (k) (blue), S2{k) (orange, dashed) are plotted in Fig.[5]for 7 = 4.52, po = 1, K = 
1.81 as functions of k/kp. We can analyze the limiting behavior as k — > of these 
prefactors using their forms 

simrfl L 9(5uj) + simrfl R 6(-Suj) 2nSi(k)6uji iR +i lL - 1 



S(k,w) 



S{k,uj) = 9{5u) 



-Slu) 

shi7r(/iL + fj R ) r(/I R + fL L )(v + VdY L \v - Vd^ 11 

\Suj\ = \uj-ei(k)\ < ei(fe), 
27rS 2 (k)8uj^ + f :iL - 1 



r(A« + h)(v + v d )^ \v - v d \^ ' 

and the universal relation for S(k, to) given in Ref. 

m*K „, k 2 



5lu = uj - e 2 {k) < e 2 (k), 



S(k,w) 



v 2m* 



\lo — v\k\ 



(199) 



The exponents P^r(L) = (-^(i^M + V K/2 + l/2\/K) and m* is the effective mass. 
In the limit that k — ¥ we have fin — ¥ 1, /ij, — ► 0. Then from the above relations we 
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obtain 

YvcaS 2 {k) = S l {k) = ^-. 

k— >0 Z7T 

This asymptote is indicated by the red dotted line. 
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Figure 4. (Color online) Results for the Lieb-Liniger model of ID bosons: 
2n 2 Ai (dashed black), 5127r 10 A 2 (solid green) and B (dot-dashed blue), — 8-7r 2 Bi 
(dotted red) as functions of the Luttingcr liquid parameter K. In the limit of 
strong interaction (K — > 1) our expressions for Bo and B\ agree with the known 
values [7], while Ai — > 1/2-7T 2 , A2 — > are in accordance with the density correlator 
of the free Fermi gas. We also match Bo in the weakly interacting regime (K ^s> 1) 
to Popov's result (dashed line) [6], and show some numerical results (crosses) |10| . 



Numerical results for A ± (k), A ± (k) as functions of k are obtained from 
Eqs. dH5j), ([Ml) and (figg]) . and are plotted in Fig. ©for K = 1.81. 



■52 
(200) 
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Figure 5. Plot of 27rS'i (blue), 2ttS2 (orange, dashed) as functions of k/kp for 
K = 1.81. As k/kp — > they approach K indicated by a dashed red line. 



5. Prefactors of equal time spin correlators of XXZ model with finite 
magnetization 

5.1. Introduction 

Let us now turn to the XXZ magnet, whose Hamiltonian is given by equation ([5]). 
For the sake of brevity, our discussion here will follow that of the Lieb-Liniger model, 
and omit extensive discussion of the various steps, which are similar. 

All eigenstates of the XXZ magnet can be obtained from the Bethe Ansatz 
|40j . Assuming periodic boundary conditions, the quasimomenta k appearing in the 
wavefunction are parametrized in terms of rapidities according to the definition 



sinh(A + iC/2) 



2 at an 



tanh A 



(201) 



sinh(A-iC/2)' " tanC/2 

in which the anisotropy A (we restrict ourselves here to < A < 1) is parametrized 
by C = acosA with £ G]0, 7r/2[. The Bethe equations for an eigenstate with M down 
spins are 

sinh(A :) - + i(/2) N -^4 sinh(Aj — A^ + i() 



sinh(A J - </2) 



n 



sinh(Aj — Afe — i£) ' 



.? = 1,...,M.(202) 
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or more conveniently in logarithmic form, 

M 



6=1 



2 *N 



(203) 



in which I a is half-odd integer for N — M even, and integer for N — M odd, and the 
kernels are defined as 

0„(A) = 2 atan ta f\* ■ (204) 



tan(nC/2) ' 



Such a state has energy given by 

M 

E = JJ2 



sin ( 



o=i 



cosh2A a — cosC' 



with momentum 



M 



sinh(A a + i(/2) 



N Z~i a 



(205) 



mod 2tt. (206) 



_sinh(A a -zC/2)_ 

In contrast to the Lieb-Liniger model, there exist both real and complex-valued 
solutions to the Bethe equations for the XXZ model (we refer the reader for example 
to [67, 58 for details). For simplicity here, we will confine our discussion to purely real 
solutions only, similarly to the discussion contained in |61j (our equations are however 
directly applicable to eigenstates with negative-parity one-strings; higher strings need 
to be treated following the logic of [68]). In essence, we consider a finite field leading to 
a generic incommensurate filling, for which the ground state is given by a filled Fermi 
interval of M < N/2 adjacent quantum numbers/rapidities centered on the origin. 
Starting from such a state, excitations can be constructed using a similar logic as for 
the Lieb-Liniger case, namely by removing/ adding quantum numbers and applying 
particle-hole excitations. 

To obtain correlation prefactors, our starting point is once again the finite- 
size representation of matrix elements of the considered spin operators. For the 
Heisenberg chain, these matrix elements were obtained within the Algebraic Bethe 
Ansatz framework in [48) . We thus start here by writing down in our conventions all 
the formulas that we are going to need in further calculations. The form factors for 
the operators S z , S* can be written as 

smh(fXj + «C/2) 



i(M M i^i{A} M )i=^^p A -p (1 n 



x Y[ I sinh 1 (fj lj - [i k - 
j>k 



sinh(Aj +</2) 
- iQ sinh - (Aj 

2P ab (M,{A})) 



A, 



i0\ 



detl/ 2 d>(M)detl/ 2 $({A}) 



(207) 



\({n} M+1 \s-\{\} M )\ = vm. 



,1/2 



(0 



n 



M+l 



sinh(/ij — «C/2)| 



I\T=i I sinh(Aj 



</2)| 



M+l 



A/ 



n i smh- 1 ^- - - *'oi n i sinh_i (^ - - oi 



i>k 

Idet 



j>k 



M+l 



det^ 2 +1 <i>(M)det^$( { A}) 



1/2, 



(208) 
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In the above expression {A} and {fi} fulfil the Bethe equations (|202[) . 

The norm of the eigenstates is given by the determinant of the Gaudin matrix 

[MUSSES], 



* ab ({A}) = 5 ab 



N- 



sinC 



E 



sin2C 



sinh A a + sin (/2 sinh (A a — Afc) + sin ( 



+(i-M- 



sin2C 



-5-. (209) 
sinh (A a — Aft) + sin £ 

The matrices H* b ({h}m, {X}m), H~ b ({fi} M+1 , {X} M ), P a b({{J>} M , {A}m) are 
given by the following equations 

TT M • \,1 \ -f\ sinh JV (At,+iC/2) t-tM . , , . .*n 

uz , , , r , , , ll j¥o smh( Mj - Afc - »C) ~ 8inh * At - iC / 2 U^a wahW - X b + <) 

smh(/i a - Afc) 

Hab (Mm+1,Wm) = i^afc ({/*}**, {X} Al) , b=l,...,M, 
H a,M+l ({M}M+1,{A}M) 



sinh /i a + sin 2 £/2 



]lm=i sinh(A m - At - iQ 



sinh 2 \x a + sin 2 £/2 

This form of the matrices is not really appropriate for taking the thermodynamic limit. 
As it is shown in Appendix E| their determinants can be recast into the following form 
given by Eq. (|E.33Jl . 

M 

dot (H* b - 2P ab ) = (-1) M sinh- M (zC) I] sinh (^ - X i - <) d ?H S ab + G z ab ) 



fr( 1+ sinh^Aj + </2) jT^ sinh(^ m - A, + 

fj ^ sinh^A, - iC/2) n„ =1 sinh( Alm - A, - iQ J 



¥r sinh(A m + </2) rTj>i sinh(^- - m) sinh(A J - A,) 
x 11 



L sinh( Mm + iC/2) J]** sinh( A i l - A, ) 

where the M x M matrix G z ab is given by Eq. (|E.32j) 



(210) 



sinh«) Ilm=i sin h(Afc - /x m ) ■^j sinh(A TO - Xb - iQ 



A I 



Uab ~ , „-/-/o\7^A7 _;_,_/ x " 11 



sinh(Afc + iC/2) H^ =1 ,fc sinh(Afc - A m ) ^Jj sinh(/i m - A 6 - <) 



A/ 



X 1 + 



siuh (Afc + </2) yr sinh(/i m - A b + iQ 
sinh Ar (Afc - iC/2) ^ sinh(/i m - A b - <)^ 
sinh(Afc + 3/2iC) sinh(A b - if/2) 



A/ 



m(A a - Ab - <) sinh(A a - A fc + iQ 
sinh(A m — «C/2) sinh(i£) \ 



2 TT plllll l^m ~ L ^l*l ainiiy^; 1 f211') 

JU sinh( Mm - iC/2) sinh(A a - »C/2) j ' 1 J 

Similar calculations are also possible for the second determinant and lead to Eq. (|E.42|) 

Af+l Af 

det H- b m, {A}) = [] [] - A, - <) det (5 ab + G£) 
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sinh Ar (Ai + «C/2) "j-j 1 s hih(jU m — K + iQ | 



n 



sinh Ar (A l - </2) sinh(/i m - A, - iQ 



I sinh(Aj + </2)| 2 n -m 1 sinh( Mj - M j) sinh (A, - A 



M 



nfIt 1 |sinh( Ml+l C/2)P 



IIS 1 Iltei sinh(/ii - A,-) 



where the M x M matrix G b is given by Eq. (|E.41[) 
1 



Hm=i sinh(A b - fi m ) n*Li sin h(A m - A fc - iQ 

sinh(A b + zC/2)| 2 n^ =1#6 sinh(A h - A m ) sinh(/x m - A 6 - <) 

x ( 1 sinh jv (A fc + zC/2) ^yi 1 sinh(/i m - A b + <) \ 
X y sinh Ar (A 6 - »C/2) sinh(/i m - A 6 - <) J 

x / sinh(A b - iC/g) sinh(A b - 3</2) g \ ^ 
V sinh(A a - Aft - iQ ' J ' 

In this new form the expressions for the form factors (Eq. (|207p and Eq. (|208[1 ) read 
w-^s \det M {Sab + G z ab )\ 



M 



detif $({ M })detl/ 2 <f>({A}) 
sinh(^tj — Aj — «£) 



sinh 1 ^ 2 (/i J ; 



»C) sinh 1/2 (A l - Xj-iO 



n,>j sinh( ( u J - jiij) sinh(A i - Aj) 



sinh(^ - Aj 



><n 

i=l 

\m M+ i\s-\{x} M )\ 



sinh iv (A l + </2) UZ=i sinh( Mm - Aj + iQ 



sinh (Aj - </2) Il m= i sinh(/i m - A; - <) 

N\ sinHK) \6 q , Pft Px J de M^ + ^)l 

A deti/ 2 +1 *(M)deti( 2 $({A}) 



(213) 



Il^t 1 UjLi sinh(^i 



Xj - »C) 



n- - +1 sinh 1 / 2 ^ - Mi - iC) Il5 Binh 1/2 (A, - A, - <) 



nf£i|sinh(A i+ zC/2)| 
n^lsinhO^C^)! 

n 



Il^t 1 sinb-Oj - fii) UjU sinh(A i 



A,) 



rii=i n ? =i smh(^ 



sinh" (A, + «C/2) ^tt- 1 sinh(/x m — Aj + iQ 



n 



sinh (Aj - </2) ^ sinh(/i m - Aj - iQ 



(214) 



5.£. Thermodynamics 

In the continuum limit, similarly to the Lieb-Liniger case, we can define the root 
density function p(X) and backflow shift function F by [33] 

p(A) = oi(A) - r dw(n)a 2 (X - fi), (215) 

J-q 

F(X\^,„-) = -f d^W\ ^{X -„)- MX ~ " +) ;_ MX ~ (216) 



2tt 



PexcMVi) = PN{Xi) + -r= i^(Aj) - F N (Xj)^-^ , (220) 
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The kernels appearing in these equations are given by (|204l) and 

2tt dX 7r cosh(2A) — cos(nC) 

The rapidity shift function is related to the order 1/N change in rapidity positions 
under the addition of excitations by the definition 

»<= x >--mj + ° {1/N2) - (218) 

Looking carefully at the finite size corrections one obtains in the same way as for the 
Bose gas (compare with section l4~2l) the following equations 

MAO = oi(A<) - f dfip N (fi)a 2 (Xi - At) + G(l/N 2 ), (219) 

r ,\i -i f q a tp t i + (\ \ ^(A-Mo) -02(A-Mo) 
-F/v(A|/i , /j ) = -/ diiF N (fi\^,fi )a 2 (A-/i) ^— — 

5.5. Thermodynamic limit 

5.3.1. Gaudin determinant The Gaudin matrix for the ground state given by the set 
of rapidities {A} reads 

*ab ({A} M ) = S ab ^2vr^ ai (A a ) - 27r^a 2 (A - X k ^J + 2™ 2 (A a - X b ). (222) 

Its determinant can be written in the thermodynamic limit as 

M 

det$ ab ({A} Af ) = [](27r^ jOA r(A. i ))det(l + a 2 ) x (1 + 0(1/N)) . (223) 

i=l 



--T7 a 2 (A - p, ) I — — r- -TTT- - a 2 (A - Mo i 



By analogy for an excited state given by a set of rapidities {/i} we get 

M 

det$ ofe ({/i} M ) =Y[(2ttN P 

exc,N (//*)) det (1 + ct 2 ) 

i=l 

= ft exp (/' dX (™ - FnW ^) ) IS* $Qb ({a}m) • (224) 

5.3.2. Pref actor Let us start by reorganising the prefactor of the S z form factor 
Eq. (|213l) . We consider first the thermodynamic limit of the following expression that 
resembles the Bethe equations Eq. (|202[) 

Mo(AQ = Smh !! (A * + K,2) n r 1 Sinh(M? " = ^ + K) ■ (225) 
sinh Ar (A l - K/2) n*f =1 sinh^ - A, - iQ 
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By separating the excited rapidities, using the Bethe equations and taking the 
thermodynamic limit of the smooth part one obtains 

( f q 27ri f q F 2 (u) \ 

M (Xi) = - cxp [2iri J d^F N (n)a 2 {n — \) H — — J c ^^ a' 2 (> - A 4 )J 

x yr sinh^-A. + zQ sinh( Mm -A t -»C) + 0{1/N 2 } _ (226) 
fj^ smhiutt - Xi - i() sinh(/i m - A, + i() 

This can be recast into a more illuminating form using the following relation 

sinh(/i ± - Xi + i() , ± . 

— — — — = cxp (uf) 2 A, - ^)) 

= exp (i^ Mo*) + i^{Xi rf) * N( fl ) + 0(1/N*), (227) 

and together with the equation for F(X) we obtain 

MoiX,) = - cxp {-2mF{Xi)) 

( 2m F N (Xj) ^ , + _ > 27ri J^(Aj) , \ 

X ° XP ( AT MX~) 2^ M A * - ^) - a2 ( A « - ^)) + 2iV^) ° l( J 

^G*^ jy^-^ (^^M-n^)^)) . (228) 

We can make an even further simplification if we consider equations for p'(X) and 
F' (A), which read 

P '(X) = a[ (A) - f d W (»)a' 2 (X - ft), (229) 

■'-9 

F'(A|/x+ /O = - / dnF(p\n+, /i-)a' 2 (A - /i) - a 2 (A - /i+) + a 2 (A - /O- (230) 



9 



and finally 
M (A,) = - cxp 
When applied in the prefactor one obtains 



(231) 



" / sinh jv (A t+t C/2)n^ 1 Binh(M m -A a + <) \ _ " 

S I S in^(A 4 -</2)n^ 1 sinh (Alm -A 4 -<)J " ^ + M ^ ^ 
= f[ (2e** sin^F(A 2 )) cxp (tt ^ dXF(X) cot (ttF(A)) (V(A) - ^(A)^ 



In cf> we gathered all purely phase factors. As we are interested only in the norm 
squared of the form factor, such a phase does not contribute to the final answer, and 
we discard it. Llet us define the following notation 

Mr = fl —TB Slnh( ^ = Aj 'J? , (233) 

ijii sinh 1/2 ( Mi - m - i() sinh 1/2 (A, - A, - i() 



M 2 = TT sinhl/2 (^-^)sinh 1 / 2 (A J -A a ) 
• f 1 , sinh(/ii-Aj) 
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and take the thermodynamic limit of the diagonal part of the denominator that is 
missing in M 2 as compared with the expression for the form factor, which reads 

n 7 ■ x n tt sinh (Pi _ X i) TT N PN(^i) 

^smh fc-A^n^^n-— 

TT F N (X l ) yr NpyjXi) 

t\ N PN (Xi) sinh(/i+ - Xi) l\ F N (Xi) ■ [ ' 

5-4- Partial results 

5.4-1- Sq form factor Let us now summarise the results of the last section. Taking the 
thermodynamic limit of some of the expressions we obtain using Eq. (|213[) . Eq. (I223[) , 
Eq. (|224l) and the above discussion about the prefactor 



\(Wm\S \{X}m)\ = -z-Sq^-p* 11 - 



i=l 



Nsmh(fil - Xi) J 

P'ip) 
Pip) 



cxp ^ alpF(p)cot(nF(p) (f'(h) - ±F(ji) 
exp(-i£^(F'W-F(M)^) 



Calculation of the thermodynamic limit of the last line of this expression is contained 
in |Appendix F| 

5-4-2. S~ form factor In the case of S~ Eq. (|214p similar manipulations are possible 
and their result are shown below. Before that however we need to take a look at the 
shift function and labelling of rapidities. We start with the shift function. In the case 
of creation/annihilation of a rapidity it takes the following form (by analogy to the 
Bose gas, see discussion above and below Eq. (|89|) . Eq. (|90|) and |21| ) 

F + (A) = F (A; {/i+}„+i, {pi) n ) - Z(X)/2, (237) 

F_ (A) = F (A; {p+}n-i, {^}n) + Z(X)/2, (238) 

where Z(X) is a dressed charge given by the solution to the following integral equation 

m 

Z(X)+ I dpa 2 {X~ p)Z{p) = 1. (239) 

J-q 

Note that the convention here is a bit misleading since acting with the operator S~ 
creates an extra excitation as compared with the reference state and therefore the 
shift of the ground state rapidities is given by F+. The opposite holds for the 5+. 

There is therefore always, by definition, one excited rapidity, which we will call 
A*n+iJ where n is number of additional excitations. In order to keep the notation 
simpler we relabel the set {pi]M+i as follows. We single out the excited rapidity 
p-n+i an d relabel the rest as {/2i}jwf 

{Hi}m+i = {^}m U (240) 
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Let us repeat that such a decomposition is always possible for the excited state in the 
S~ form factor. 

To show how it works in practice let us consider in details the following expression 
from the prefactor of the S~ form factor 

1/2 



(nf^ 1 sinhOj - Mi) smh(Xj ~ A,)) 



T, 



extra 



x Mo 



UZi 1 Ilfei sinh(/ii - Xj) 



lj=i vr. • -j , Il l JV Ii sinh(/i i - Xi 

where we defined a similar expression to M2 Eq. (|234[) and T extra 



m 2 =n 

Textra 



sinh 1 — fii) sinh 1/ ' 2 (A. ) — A;) 



sinh(/i,; — Aj) 



M 

n 



sinh (p n +i - /ii 



Ai- sinh(^ n+ i - Xi) 
The product in the denominator can be rewritten as 

M M AT l\ \ n 

nsmh-^-A^n^n 



L\ t\ Np(Xr ) s inh(rf Xr) 



(241) 

(242) 
(243) 

(244) 



Actually, terms M2 and M% are equivalent, the only difference being that we have to 
use the modified shift function F + (X) in the latter. In the same way we obtain the 
following equality 



IlfiV UjLi sinh (^ - A 3 - »C) 



n- T sinh 1 / 2 ^ - N - iC) n- ■ sinhV^A, - X j iQ 



Mi 



sinh(— iQ 



1/2 M 
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(245) 



Keeping in mind the relabelling of the rapidities and following the same logic as for 
S* we can obtain the following partial result for S~ given in Eq. (|214j) 
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(246) 



Again the last line is the only line where we still need to take the thermodynamic 
limit. As mentioned before calculations of M\ t 2, M\ t 2 are similar to each other. The 
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difference is hidden in the shift function. In the case of S* operator the excited state 



rapidities are connected to the ground state by 
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while for S„ 



flj — Xj 
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N P {X t 
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(248) 



Therefore we will perform calculations only for S* and simply adapt the results by 
changing the shift function in the final formula for S~ . 



5.5. Results 

5.5.1. fo rm factor Collecting all the partial results from previous sections and 
|Appendix F| (Eq. (I2lt6")) . Eq. (|F.12j) . Eq. (|F.47[1 ) we obtain a final expression for the 
form factor of the operator S* 
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where we used the following property of the Barnes functions 

G(l + F(q))G(l - F(q))G(l + F(-q))G(l - F(-q)) x 

dXTrF{X)F , (X) cot( 7 rF(A)) N ) = G 2 (l + F(-q))G 2 (l - F(q)) (2ir) F{q) - F{ - q) 
1 J 



x exp 



Umklapp form factor If the excited state is formed by m consecutive umklapps, then 
we have a number of simplifications. First of all we do not make any mistake if we set 
(i = 1, . . . , m) 

Pi = -q, (250) 
pf = q. (251) 
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Then going term by term we have 
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As calculated before the Thole term takes the following form 

xexp^P./ ^J^J , ,254, 
and Tparude is given by 

w - «,»,(,»•■-» ft x «p (_ m p + /_; rf ,_j^L s ) ( . 2 55, 

The T cross term reads 

T cross = (iVp(g) sinh(2 (Z ))"" l2+m G 2 (m + 1). (256) 
Combining together Eq. (|254l) . Eq. ([255)l . Eq. (|256| we obtain 
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G(m + l-F(q))G(m + l-F(-q)) 

T m (F(q))T m (l- F(q)) ' 1 ' 

For the umklapp excitation there is again a relation between the shift function F{\) 
and the Luttinger liquid parameter K. Function F(X) fulfils an analogous equation as 
in Eq. (fT78|) . 

F(X) = — T dp (F(p) — 1) a 2 (X — p). (258) 

J-q 

This equation is solved by F(X) = 1 — Z(X), where Z(X) is a dressed charge defined 
in Eq. (|239|) . The value of the fractional charge at the Fermi boundary is connected 
again with the Luttinger liquid parameter 



Z(q) = VK; (259) 
and therefore for the m-th order umklapp we obtain 

F(X) = —m (Z(X) — 1)) . (260) 
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where we used the following relation 

1 = sinh(2a)-zsin(2b) 

tanh(a + if)) cosh(2a) — cos(26) ' 

and the relation between the form factor of the umklapp state and the prefactor of 
the correlation function (Eq. ©, Eq. ([TT])) 
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5.5. 2. S q form factor Gathering results from previous sections Eq. (|246p . Eq. (|F.12p . 
Eq. (|F.47I) we obtain 
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x f[ f ^ffi 1 "*!"!^ G(l + F + (-q))G(l - F + (q)). (264) 

Umklapp excitation Let us start with calculation of T extra defined in Eq. (|242[) . it 
gives 

T t _ tt sinh (/Xn+i - Ai) = yr sinh (m„+i - -A- sinh(/i„ + i - ft+) 
extra ii sinh(Mn+i - A,) Al sinh (/x„+i - A,) Al sinh(// n+ i - ft) ' 

The first product resembles T part icie and indeed for the umklapp excitation it reads 
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where we put the rapidity Mn+i exactly at the right Fermi point. The second product 
gives 
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where we used that (see Eq. (fTMl) . Eq. ((237)) ) 

,F + (A) = mF(A|Umklapp) + F(A|g) - Z{\)/2. (270) 

and 

F(±g|«) - Z(±q)/2 = I ± 1 (l - J=Y (271) 



Using the relation between umklapp form-factor and the prefactor of the correlation 
function (see Eq. ©) 
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we obtain 
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6. Conclusion 

We utilized a general method for obtaining non-universal prefactors of correlation 
functions from field theoretical considerations. The technique relies on the existence 
of a well defined relationship between lowest energy form factors of operators and 
the associated prefactors of their correlation functions, see e.g. Eqs. (ThO - fTS)) . For 
example, in the case of equal time correlators, such a relationship can be demonstrated 
by an effective field theoretical description of the system as a Luttinger liquid, when 
finite system size is properly accounted for. In fact such relations were already known 
in a few cases [64] . Moreover using the three subband model of a mobile impurity 
in a Luttinger liquid, we also obtained expressions for prefactors of dynamic response 
functions involving lowest energy form factors, see Eqs. (|2T|) . (j3TJ)) and ([53")) . 

The universality of the field theoretical description of interacting one-dimensional 
quantum systems of bosons, fermions and spins, allows us to apply the relationship 
between the correlation prefactors and form factors to a wide variety of systems. 
It should be emphasized that the method is quite general, and does not rely 
on integrability or on any other special features of a particular system under 
consideration. 

In this article, we focused on three integrable models: the Calogero-Sutherland, 
Lieb-Liniger and XXZ models. We utilized the connection between prefactors and 
form factors in the case of the aforementioned systems, where explicit expressions for 
finite-size form factors were already available. We demonstrated how to properly take 
the thermodynamic limit of such form factors, which contain non-trivial power laws 
in system size, to obtain numerically tractable analytic expressions for non-universal 
prefactors of their correlation functions. 

We note here that some results for prefactors we obtained do not represent new 
results. On the other hand, it was our aim to demonstrate a more direct way to obtain 
them, with a wider range of applicability, e.g. we obtained the prcfactor of the first 
oscillatory component of the equal-time density-density correlator for the Calogero- 
Sutherland model without employing the replica method or multiple integrals |12j . 
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similarly our results for the correlation prefactors of the XXZ model did not require 
use of a master equation, or detailed asymptotic analysis of determinants [13] . 

We summarize the results for prefactors of correlation functions in the following 
sub-sections. 



6.1. Prefactors of equal-time correlators 



The equal time correlators of ID quantum liquids can be expressed as an asymptotic 
series with unknown prefactors, see Eqs. (JlJ - (J2J), ([S]) and ©. From the analysis 
of the finite size scaling of the field theoretical description of ID quantum liquids 
one obtains a correspondence between the prefactors and form factors of the density, 
creation, annihilation and spin operators, see Eqs. ([T4]) - (IT8|) . We present below the 
prefactors of the equal time density and field correlators of the bosonic Lieb-Liniger 
model followed by prefactors of the spin-spin correlators of the XXZ model. 

For the first prefactor A x of the oscillatory terms in the density-density correlator 
(see Eq. (JTJl) of the Calogero-Sutherland model we obtained 
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in agreement with the results of Ref . [T2] which were obtained using the replica method. 

For the prefactors of the oscillatory terms in the density-density correlator (see 
Eq. ([!])) of the Lieb-Liniger model we have from we have from Eq. (TH)]) . (|170[) and (|181|) 
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with F(X) = m{\ — 2irp{X)), where the quasiparticle distribution function p(A) is 
given by Eq. (|60[) , q is the edge of the quasimomentum distribution and satisfies 
2np(q) = K, the Luttinger parameter, po is the density and c is the strength of the 
interaction between the bosons, with the associated dimensionless parameter 7 = c/po- 
G(x) is the Barnes G function defined in Appendix A, and the symbol P± defined in 
Appendix B is an instruction to evaluate the integral in the principal value sense 
where the singularity occurs at the edge of the range of integration. G is obtained 
from Sec. 4.5, Eq. (|163l) while K(X) = -^r^i- The two determinants are meant to be 
evaluated as Fredholm determinants, see Appendix D and Ref. |75| . 

For the prefactors of the oscillatory terms in the Green's function (see Eq. (j2J) of 
the Lieb-Liniger model we have from Eq. ([T4|) . (|186j) and (|191j) 
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The operator G z is given by Eq. (|F.56[) and a 2 (A) 



i(2C) 
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A 



where F_(A) = F m (A|Umklapp) — F(X\q) — Trp(X) is the modified shift function defined 
in Eq. (19TJ1) . if is the Luttinger parameter, q the edge of the distribution of ground 
state quasimomenta, po is the density and c is the strength of the interaction between 
the bosons, with the associated dimensionless parameter 7 = cj pa. G(x) is the Barnes 
G function defined in Appendix A, and the symbol P± defined in Appendix B is an 
instruction to evaluate the integral in the principal value sense where the singularity 
occurs at the edge of the range of integration. H~ is obtained from Sec. 4.5, Eq. (I169[) , 
while K(X) — ^ 2c ^ . The two determinants are meant to be evaluated as Fredholm 



determinants, see Appendix D and Ref. [75] . 

For the prefactors of the correlations functions of the xxz spin chain we obtained 
analogous expressions. The prefactor of the S Z S Z correlation function (see Eq. (|8]l) is 
given by (cf. Eq. (!?]), Eq. pHTjl ) 



(277) 



,, cosh(2(g - A) - cos(2C) , 

where N is the number of spins, q the edge of the distribution of ground state 
quasimomenta and p(q) is the density at the edge of the distribution. The Luttinger 
parameter K is given by K — Z 2 (q) where Z(X) is the dressed charge defined in 
Eq. (|239[) . The interaction parameter £ is connected with the anisotropy A by 
A = cos(C). The shift function is given by F(X) = —m(Z(X) — 1) and G stands for 
the Barnes function defined in Appendix A. The symbol P± is defined in Appendix B. 



The two 
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determinants are meant to be evaluated as Fredholm determinants, see Appendix D 
and Ref. [75]. 

For the prefactor of the S + S~ correlation function (see Eq. ©) we have (cf. 
Eq. (USD, Eq. (gTJD ) 
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(278) 

The meaning of all the symbols in the above formula is the same as in the formula for 
D m , with the exception that the shift function is given by f+(A) = rni 7 '(A|Umklapp) + 
F(X\q)-Z(X)/2. 

The numerical evaluation of Eq. (|277|) . Eq. (|278[) shows that only prefactors with 
small value of m (m = 1,2 for D m and m = for _E m ) have significance. These 
prefactors are plotted in Fig. 6-9 as functions of the anisotropy parameter A and the 
filling M/N. 

We note that similar results for equal-time correlation prefactors were obtained 
in Ref. [14) . but establishing their equivalence to ours proves non-trivial due to the 
appearance of multiple integrals. 

By applying the techniques described in Sec. 2 above to the prefactors of 
singularities in dynamic response functions like the spectral function and density 
structure factor, we obtained a correspondence between form factors of the creation 
and annihilation operators and the prefactors of the response functions. We present 
results for the prefactors of the dynamic response functions of the Lieb-Liniger model 
below. 



6.2. Prefactors of the singularities of the density structure factor 

The density structure factor S(k,ui) defined in Eq. (|2"TT) shows singular behavior for 
lj s» \s(k) = A(fc 2 — k%)\, with prefactor S'csm(^) when < k < 2kF- The prefactor is 
given by 

In the case of the Lieb-Liniger model, the density structure factor S(k,ui) also 
shows singular behavior for u £\{2){k), i.e. near Lieb's modes [21] [38], and the 
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behavior there is described by Eqs. (j2"5")) and (j2"6")l . For the prefactors S\(k) and S%(k) 
in those equations we obtain from Eqs. (1271) and (11701) 
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Note that the momentum dependence in the above expression is contained in the shift 
function F(X) = F(X\/j, + ) — F(X\q), defined in Eq. (fBTj) as well as the terms directly 
involving the quasimomentum corresponding to particle excitation (see Eq. I|49p). 
We may relate the physical momentum k to the quasimomentum p + using [23J [5T] : 

k = fi + - irp Q - f 1 d\9{n + - A)p(A), (281) 



where, 8(X) = i log ^ffrj^ , and q is the quasimomentum at the edge of the 
distribution; po is the density and c is the strength of the interaction between the 
bosons, with the associated dimensionless parameter 7 = c/po. G(x) is the Barnes 
G function defined in Appendix A, and the symbol P± defined in Appendix B is an 
instruction to evaluate the integral in the principal value sense where the singularity 
occurs at the edge of the range of integration. G is obtained from Sec. 4.5, Eq. (|163j) . 
while K{X) = c ^ X i ■ The two determinants are meant to be evaluated as Fredholm 
determinants, see Appendix D and Ref. [75] . 

Similarly, we get the following expression for ^(fc): 



S 2 (k) 



c 2 VK 



-(F 2 (q)+F 2 (-q)+2F(q)+2) 



exp 



exp 



exp 



9 



2F(X)(q-p,- 



dX- 



dp 



ic)(X — q + ic) 
F(X)F(p) 



(27T) 



/„„ /, „, / G{l + F{-q))G{l-F{q))T{\~F(q)) Y 



r(i + F(/x-))r(i-F(/x-)) 

F 2 (^) 



(A - v + ic) 2 
1 ri 

x — p~ /q 2 



dx 



2F 2 {qx) 



dX I d/i 



i x 2 - 1 

F(X) - F( P ) 
X — p 



- R 



1 



x-l 



Det 2 (l + G) 



Det 2 1 



K_ 

2tt 



(282) 



Here too the dependence on momentum k is carried by F(X) = F(X\q) — F(X\p,~) and 
\i~ , the quasimomentum of the hole excitation (see Eq. (14^1) ). This time we relate the 
quasi momentum of the hole, pT , to the physical momentum, k using 

- k = pT - npo - [ dX6(p- - A)p(A). (283) 
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6.3. Prefactors of the singularities of the spectral function 

The spectral function A(k, lS) defined in Eq. ([2"2")l also displays singular behavior in the 
vicinity of Lieb's collective modes [55] • This behavior is described by Eqs. (J25J), (pJTl) 
and the text in that section. There are four associated prefactors defined in Eqs. ([28l) 
and (|3ip for which we obtain analytic expression below. 



6.3.1. A±(k) We obtain A + (k) from the creation operator form factor for a Bosc 
gas with a high energy particle, and obtain A- (k) from the annihilation operator form 
factor for a Bose gas with one high energy hole. These functions give the prefactors 
for the spectral function singularities at ei(fc) and — £2(fc) respectively. We must first 
relate the quasi momentum of the particle to the momentum k of the excited bosonic 
state as follows [53J [5T] : 

k = p+ - 7T Pa - r dX9(p + - X)p(X), (284) 

J-q 

and for the quasi momentum of the hole 

- k = pT - Trpo - r d\0(fj,- - X)p(X), (285) 

J-q 

where 9(X) = i log f ffzj ^ , and q is the quasimomentum at the edge of the distribution. 

Particle-hole pairs in Eqs. and (IM1) are defined with respect to ground states 
of TV ± 1 particles, so to obtain a state with a single high energy particle (hole), a 
particle-hole pair needs to contain a hole (particle) at q. Such a procedure results in 
the following expressions for shift functions F± (A) 

F+(Xj = F(X\p + ) + np(X), 
F_(X)=-F(X\ l x-)-7rp(X). 
From Eq. ([30]), ([PM]) and (fl"87]) we have 
1 



(286) 



A±{k) = 



exp \ P± 
exp 



q^K 



(27T) 



F ± (q)-F ± (-q)-2 



(G(l + F ± (- q ))G(l-F ± (q))Y 



dp 



x- 



1 



dx 
-1 x 



2F ± (qx) 



F±(A)-F±0i) 



A 



rfA 



exp 



±2 



dX- 



F±(X) F ± {p) 
(A - p + ic) 2 

F^(X)(p ± - q) 



Det 2 (l+H-) ((q-p 



±\2 



Det z 1 - 



K 



(A — p ± + ic)(X — q + ic) 



(287) 



6.3.2. A±(k) Similarly, we obtain the prefactors of the spectral function singularities, 
A + (k) and A_(fe), near the Lieb modes e%(k) and — e\(k) from the form factors of the 
creation and annihilation operator respectively. We can obtain the prefactor A+(k) 
from the creation operator form factor of a system with a high momentum hole and 
two particles at the right quasi- Fermi point. From the annihilation operator form 
factor of a system with a high momentum particle and two holes at the right quasi- 
Fermi point, we obtain the prefactor A_(fe). We again relate the quasimomenta of 
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the particle and hole to the momentum k using, Eqs. (|192[) and (|193[) and determine 
F±(A) using Eqs. (55 ) and (pi )) . 

Particle-hole pairs in Eqs. (j59")l and are defined with respect to ground states 
of N ± 1 particles, so to obtain a state with a single high energy particle (hole), a 
particle-hole pair needs to contain a hole (particle) at q. Such a procedure results in 
the following expressions for shift functions F±(A): 



F+(X) = 2F(X\q) - F(A|/0 + ttp(A), 
F-(A) = F(A| M + ) - 2F(A|g) - 7rp(A). 
Then using Eqs. (| 186[) and (|187[) we get the following expressions: 



(288) 



- (F* (q)+Fl (-g) -4F + (g)+4) 



p( M -)( M -/9-l) 4 
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(27T) 
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i 2 P + (gx) 
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-i /o 



F + 

1 

2 



i 4F+(gx) 
dx ; - + 2 

-l 



x-1 



F+{q)) \ Dct 2 (l + g+) 
/ ^(l-f) 
F + (\)(q-ii-) 



dX- 



(A — + ic) (A — o + ic) 



dX 



dp, 



'P + (A)-P + (/i) > 



<-/// 



« F+(A) F+Qu) 
uA ~ 

(A-^ + ic)^ 
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p( M +)( M +/ -l)4 



gVtf 

+ ^2 



-(Ff (g)+F=(-«)+4F_(«)+4) 



(2tt 



F_(g)-F_(-g)-2 



Det 2 (l + H-) ((q - ,i+) 2 + c 2 ) / G(l + F_(-g) )G(l - F_(g) )r(l + F_(g)) r(2 + F_(g) ) 
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exp 
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I\F_(g))r(l-F_(g)) 



i 2F 2 (gi) 

dx 2 7= 
.1 X z — 1 



1 4F_ (gx) 
dx- 

-l 



x - 



1 



dA 



P-(A)-P_( M ) ' 
A — u 



exp 



exp 



-9 



, » P_(A)P_(/i) 
d/z / aA 

{X-p + ic) 2 

F_(A)Qx+-g) 
o + ic)(A — /i+ + ic) 



3 2F_ (A) 
dA- 



X 



In the expressions above, q is the quasimomentum at the edge of the distribution; 
po is the density and c is the strength of the interaction between the bosons, with 
the associated dimensionless parameter 7 = c/po. G(x) is the Barnes G function 
defined in Appendix A, and the symbol P± defined in Appendix B is an instruction 
to evaluate the integral in the principal value sense where the singularity occurs 
at the edge of the range of integration. H is obtained from Sec. 4.5, Eq. (|169[) . 
while K(X) — c 'i 2 + x i ■ The two determinants are meant to be evaluated as Fredholm 
determinants, see Appendix D and Ref. [75] . 



(290) 
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Calculations of the prefactors of singularities in dynamic response functions so far 
assumed that the field theory of mobile impurities provides an adequate description. 
The microscopic approach of this article allows one to also explicitly prove the existence 
of singularities without this assumption. If we consider the total spectral weight in 
a small interval of energy Sui in the vicinity of £1(2) (fc), we need to sum over states 
with low energy particle-hole excitations. We may use Eq. (|20|) . proven in Ref. [15] 
and the expansion C(n r (j), M-R(i)) ~ r V(Z)*'' AXm.r(z,)) for n r ^ ^ 1 to show that he 
total spectral weight scales with 5uj as 

y \C(n r ^ R )C(n lt , L )\ 

\E-e 1(2) \<8u/2 

which proves the existence of the singularity. 

To summarize, we analytically calculated "non-universal" prefactors in long- 
distance behavior of correlation functions of the exactly solvable Lieb-Liniger model 
of ID Bose gas and the XXZ model of a ID spin chain with anisotropic coupling. We 
also calculated prefactors of singularities in dynamic response functions such as the 
density structure factor of the Calogero-Sutherland model of fermions with long range 
interactions that scale as the inverse of the square of separation, as well as for the 
Lieb-Liniger model and the spectral function for the Lieb-Liniger model, and proved 
the existence of singularities within a continuum spectrum. Our results represent a 
significant step towards an analytical calculation of the full correlators. 
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Appendix A. T products and Barnes G function 

Following are the definition of the Barnes G function and it's relation to the T- function 
as well as a few special values and asymptotic expansions: 



A G(n + a+l) 

llr(» + o)= , 



G(l + a) 

G(l) = l, 

lim log(r(n)) = (n - 1/2) log(n) -n+ ^log(2n) + Ol- 



n 2 3 n 1 / 1 \ 

lim log(G(n + 1)) = — log(n) - -n 2 + - log(27r) - - log(ra) + C'(-l) + O -(1.1) 
n-Hx> 2 4 2 12 \n ) 



Appendix B. Principal Value Integrals: Interior and Edge Singularities 

ows, 



We define three distinct types of principal value integrals as follows 

r c-S 



J a X-C S ^°\Ja X-C J c+5 X-. 

J a 



C 

6 dx^- = lim ( [" K*L + f(a) \og(S)) , 
x-a s^o yj a+s x-a J 

r j d.r±- = lim ( r A ^4 - f(P) log(S)) , 



x — b s^o \J a x — b 

P± tdx f ^ b -^ = P I" W P [ b I^l (Bl) 
± J a (x-a)(x-b) + J a x-b J a x-a' 

Note that we consider all quantities being integrated to be dimensionless, and 
the bounds of integration to be real numbers with no physical units. All such special 
principal value integrals evaluated in the text are first made dimensionless by mapping 
the region of integration to the interval (—1,1). 

Appendix C. Special Property of Barnes Function 

We will use the following property of the Barnes function, 

G(l -z)= G ^ + *\ So ™ cot ^ d * (c.i) 
v ; (2tt) z ' v ; 

to simplify a group of terms occurring in the final expression for the various form 
factors. 

Now consider the term, 

exp d\iTF(X)F'(X)cot(nF(X))^ 
= exp < / dx7rxcot(7rx) > 

{JF(-q) J 

[ rF(q) rF(-q) ~\ 

= exp < / dx7rxcot(7rx) — / dx7rxcot(7rx) > 
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_ G{1 - F{q))G(l + F(-q))(2ir) F( -^ F ^ 
G(l + F(q))G(l-F(-q)) 

Thus we have 

G{\ + F{q))G(\ - F(q))G(\ + F(-q))G(l - F(-q))exp 
= G 2 (l + F(-q))G 2 {l - F(q))(2ir) F ^- F( - q K 



74 
(C.2) 

dXv:F{X)F' (\)cot(jtF(\)) 



(C.3) 



Appendix D. The Fredholm Determinant 

Given an integral equation of the form, 

<j>{x) = f(x) + X [ K(x, y)d>{y)dy, (D.l) 

we may construct a Fredholm determinant D(X) from it by replacing the integral 
above by a finite sum. The procedure to do so is as follows and is outlined in more 
detail in Ref. [75] . First we discretize the interval [a,b] into n equal parts of length, 



(b — a)/n. We may index points in the interval using Xj 



jS, and denote 



{xj),fj = f(xj),Kjk = K(xj,Xk) where j,k run from 1 to n. Using these 



definitions, the discrete version of Eq. (ID.1|) becomes, 

j = 1,2, ...,n. 



jk<Pk, 



(D.2) 



We note that the original integral equation in Eq. (jD.ll) gives rise to a system of n 
equations in n unknowns, <p\, ...,</>„. We may solve the above system by constructing 
a determinant as per Cramer's theorem of the form, [75] . [77] . 

l-XK n S -\K 12 6 ... -XK ln S 
Dt m = -W 2 i6 l-XK 22 S ... -XK 2n 6 

-XK nl 6 -XK n2 5 ... l-XK nn S 

We wish to obtain the "true" Fredholm Determinant associated with the integral 
equation in Eq. (ID.1[) as the continuum limit of Eq. (|D.3I) . One way to proceed is to 
first expand the determinant in Eq. (ID.3|) as in Ref. [77] . 



D„(A) = 1-A^ SK PUP1 
pi=i 



(-xy 



E 
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P1P2 
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1 


K 


K 


K 



and take n — > 00. This allows us to replace the finite sums by integrals. Thus we may 
define the Fredholm determinant as the entire function D(X), where 



D(X) = 1 + 



OO 

E 

n=l 



with 



dx\ 



dxo.. 



K(xi,xi) K(xi,x 2 ) 
K(x 2 ,xi) K(x 2 ,x 2 ) 



(D.3) 



K(xi,x n ) 
K(x 2 ,x n ) 



K{ ) K{x n ,x 2 ) ... K(x n ,x n ) 
The above derivation also illustrates an elementary numerical method for calculating 
Fredholm determinants associated with various integral operators. 
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Appendix E. Representation of the determinants for the xxz spin chain 

Expressions for the form factors of operators S z and S involve determinants of some 
matrices. In general the thermodynamic limit of the determinant of an arbitrary 
matrix is not a well-defined object as it may diverge when the size of matrix goes to 
infinity. The purpose of this calculations is to find different representations for the 
matrices and therefore for the determinants as well. In the final representation we 
should have only determinants of Fredholm-type. The matrices we are interested in 
are as follows 

K b {{^UA^U) = Zhh, ,M x\ 2P «o(M, {A}), a, 6 = 1, .... M, (E.l) 
smn(/z a — Ab) 



^Mg,war') = i^ { t ) . = 1 - M - 1 - (E - 2) 

1 



H; M mfi 1 ,WZi 1 )= . 7 27 ,2 . 2777^ a = l,...,M. (E.3) 

smh (fj, a ) + sin (C/2) 

In the first matrix sets {/z} and {A} consist of M elements, in the second matrix set 
{A} is smaller and consists of M — 1 elements. Matrix P a b is given by 

p rui - nf=iSinh(A m - Ab - ifl , . 

^ab({Ml, |A}J —2 . 2 , (E.4j 

smh [i a + sm C/2 

Function 3^(2, {a}) is given by 

y a (x, {a}) = [] sinhK - * - *C) - "1% /J II sinh K " * + *0- (E-5) 

The size of the products above is determined by the size of a set {a}. Whenever set 
{a.i}fL l fulfills the Bcthc equations then 

y i {a i ,{a})=Q, Vi = l,...,M. (E.6) 

The two matrices H* b , H~ b can be related to each other. Let us consider a matrix 
given by the following elements 

«.iw&.w&)- ^:':g ) , m=i,...,m. (e.7) 

Then in order to obtain i?* b we simply have to subtract matrix P a b and let {a}^ — 
Mill and {0}^! = {X}fLi- In order to obtain H~ b we take {a}^ = {/i}^ and 
{^iLi 1 = {^li^I 1 ' as weu as we se * Aw = — *C/2- Then the determinants are related 
by the following relation 



djrt2fc(M£i.{A}£r 1 ) = 

-1 



M 

= \J[ S mh( H -iC/2)\ detH^d^iX}^,-^^) . (E.8) 

All these show that we can consider first the general case of matrix H a b and 
specialise only later. The main idea is as follows. We want to extract the diverging 
part and calculate its determinant separately. 

det H ab = dct (C ab ) det (A a 6 ab + B ab ) , (E.9) 

M MM 

C ab = . r / r r - (E.10) 
smh(yU Q - Ab) 
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Then the determinants on the r.h.s can be calculated even in the thermodynamic limit. 
The first determinant is Cauchy-like, the second one can be represented as a Fredholm 
determinant. 

Let us consider all the quantities as functions of arbitrary real numbers {oa] M 
and {(3i}fL 1 . Only later we will connect them with rapidities. In order to find A a and 
B a b we simply use the inverse of C to find that 

M 

C'^Ha = A a Sab + B ab . (E.ll) 

i=l 

Such a sum can be computed by introducing an auxiliary integral of a complex 
function. This function and contour of integration should fulfill two criteria. Firstly 
the integral should be equal to 0, secondly the function should have poles that generate 
the sum. If this is the case then the sum is simply equal to the sum of all other residues 
taken with a minus sign. For the integral to be equal to it is enough for the function 
to be i7r periodic. As it will become clear later the function which would generate the 
above sum is actually odd under shift by in and therefore we introduce by hand an 
odd function f{z) which ensures that the integral is indeed equal to 0. Therefore we 
rewrite Eq. (|E.9I) as 

M 

det# o6 = TT .TV*) det (Cat) det (A a S ab + B ab ) . (E.12) 

M - LJ - M M 

and the vector A a and the matrix B a b are now given by 

M 

C^HafiOi) = AaSab + B ab - (E.13) 

z=l 

Function f(z) should fulfill one more assumption, namely it should have no poles in 
the region of integration. 



Inverse and the determinant of C ab Motivated by the fact that C a b is similar to a 
Cauchy matrix, we simply write its inverse as 

c -i _ 1 IIm=i sinh (fiq ~ gm) IIm=i sinh(a fc - p m ) , EU , 

a " sinh (^ " l£=i,*a Sinh (^ " Pm) n™=i, #b sinh(a b - a m ) ' 
One can prove that it is indeed the inverse directly from the definition. We have to 
show that 

M 

Y J C ai C a> 1 = 5 ab , (E.15) 

i=l 

which after substituting the expressions for matrix elements gives the following 
equality 

\ - -1 ]Im=i sinh(/3j - a m ) 



^ sinh(ft - a a ) sinh(ft - a b ) Y\^ =1 ^ sinh(/3 4 - /3 n 
rim=i,^6 sinh(ah — 
Um=i sinh(a b - j3 m ) 



S«b . w „ ' ■ (E.16) 
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In order to perform the summation we consider an auxiliary integral defined as 

_Ll dz Zl n^isnrh(z- am ) _ Q 

2m / 7 sinh^ - a a ) smh(z - a b ) J]^ =1 sm h(,j - p m ) 

The contour of integration 7 is chosen in such a way that the integration gives 0. On 
the other hand the integral is equal to the sum of residues. One type of residues, these 
at z — > ft gives the sum 

= ^T73 ^3 ^=1^-°") , i = l,...M. (E.18) 

z ^ a ' sinMft-a^sinMft-a^n^.sinh^-^)' 

There is one more residue when a = b then 

llm=i smh(a a - /3 m ) 
The sum of residues must be therefore 

M 

Y, Res z^ ai = -Res z ^ aa . (E.20) 
i=i 

This shows that is indeed the inverse of C a b- 

We need also the determinant of C a b- One can write 

det = TT n^i sinh (/3 J "a m )sinhK--/3 m ) ^ 



1 rim=i ^ sinh(^- - ft„) sinh(a j - a m ) « Vsinh(ft - a fe 



11^- sinhjos - ft. 



dct 



sinh(a, - ai) sinh(ft - ft ) / *l sinh^- - ft) , 

TT^ sinh(ai — ft) \ 

M UtJ ^ — J detC ab . (E.21) 

J[ j:>i sinh( aj - ai) smh(ft - ft) y 



M 



Therefore 



ntlv sinh(a,- — aA sinh(ft — ft) 
d «^= ' )/, ^ 7 P \ (E.22) 

where we used that 

detC-j, 1 - (detCafc)" 1 . (E.23) 

M V M / 

Calculation of A a and B a b Having the inverse of C ab it is rather straightforward to 
calculate A a and B a b. We follow the same logic as for the inverse of a Cauchy matrix 
and after similar calculations we obtain that 

A a ({a},{(3}) = /(ft) "— lSinh ^ TO .;/°~ <) (E.24) 

sinn(— iQ) 

sinh^ft + i£/2) n* f =1 smh(a m - ft + iQ 
AW sinh W (ft-</2) sinh( +i C) 
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and for B ab 

ni=! smh(/3 6 - a m ) l£ =ltltb smh((3 m - ft - <) 



B a6 (M, {£}) = /(ft + <) 

-/(A - iC) 



n„=i,^ sinMft - ft*) sinh(ft - ft - K) 
\C=i sinh(ft - a m ) sinh^ (ft + iC/2) 11™=!,^ sinh(ft„ - ft + <) 



' n„=i,^ sinh(ft - /3 ro ) sinh Ar (ft - »C/2) sinh(ft - ft + <) 

(E.25) 

In the formula above, we have changed the indices in the prefactor of matrix B ab such 
that it is only 6-dependent. Such an operation does not change the determinant and 
simplifies calculations. 

We have solved the problem of representing the generalised matrix in terms of 
Cauchy-like and Fredholm-like matrices. Now we specialise to the two cases of our 
interest, we start with H* b . 

Results for H* b In order to calculate H* b we have to include the extra matrix P ab , 
this is done easily because all the necessary operations are linear so they don't change 
A a nor B ab and therefore 

M 

II TT^-T d ^ ( H *bf(a a ) - 2P ab f{a a )) = dot C ab det (A a S ab + B ab + D ab ) , (E.26) 
where 

D ab = -2Y,C- i 1 P tb f(a l ). (E.27) 
Again we introduce an auxiliary integral and follow the same spirit as before to obtain 

rrM • h(\ — ) M 

{A}) = 2 n llm=lS 7^ II sinh (A m - x b - iC) 

U m =l,^a Sm H X a - A m ) AJ^ 

v/0 ,n m =i, # aSinh(A m +iC/2) ,,.^ /0 .n ro =i,^sinh(A ro -iC/2)\ 

X l /( - </2) n m = lS inh(, m + </2) - /K/2) n m = lS inh(M m -</2) J ' 

where we have already substituted the sets of rapidities {m^i} and {A^} describing 
Bcthe states for the arbitrary sets of real numbers {a^L{\ and {ft^}. If we do the 
same with expressions for A a and B ab we obtain 



]lm=l sinh (Mm - A a - K) 



A a ({n},{\}) = f(\ b ) 



sinh(— i() 

smh N (X a + iC/2) ]l™=i sinh( Mm - A a + <) ' 



(E.29) 



V sinh iv (A a - iC/2) n,„ = i sinh( M m - A a - i£) 
^(M,{A})= _ 1 Jr =1 [ " II sinh(A ro - A b - 

/(A6 + iC) /(A fc -<) 



v sinh(A a - A fc - i() sinh(A a - X b + iQ J ' (E.30) 
As the last step we transform the determinant 

M M 

[T 77-^ det (H ab f(fi a ) - 2P ab f(fi a )) = TT — - det C afc det (A a <J a6 + B a6 + £ a6 ) 
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M . 

= Iij^- ) d ^ tC ^ d ^^ + GZ ab), (E-31) 

and upon choosing f(z) = sinh(z + i(/2) we obtain 

GZ = sinh«) IIm^i sinn ( A b ~ Vm) tt sinh(A m - A b - iQ 
ab ~ sinh(A b + iC/2) \C =l ^ b sinh(A 6 - A m ) i = \ sinh( Mm - A b - iQ 

( | sinh jv (A fc + iC/2) jr sinh( Mm - A b + iQ \ 
K \ smh N (X b - iC/2) sinh( Mm - A 6 - <) J 
sinh(A b + 3/2iC) sinh(A b - i(/2) 



sinh(A a - A b - <) sinh(A a - A b + iQ 
^ sinh(A m - i(/2) sinh(»C) \ 

JU sinh( Mm - iC/2) sinh(A a - iC/2) J ' 1 - J 

Finally we can write 

M 

det^ b (M,{A}) = (-l) M sinh- M (zC) I] sinh (^ ~ A ^ ~ *C)det(<U + G* b ) 

x H / sinh^A, + zC/2) n*^i sinh(// m - A, + tQ \ 
^ sinh^A, - iC/2) nf =1 sinh( Mm - A, - »C) J 

x -pr sinh(A m + </2) TI^L sinh( Mj - - ^) sinh(A 3 - A,) 
X ii sinh (M>n + iC/2) nS- «inh(^ - Aj) 

Results for H~, Following the prescription described in the beginning we set on = fii, 
and (3i — Aj where sets {/U^Li} and {Xfti 1 } solve Bethe equations and Am = — *C/2- 
Then we obtain (a, 6 = 1, . . . , M — 1) 

{A}) = . f] sinh( Mm - Aa - iC) 

smh(— zc) - L - L 

v " m=l 

x / sinh w (A a +zC/2) p4 siDh(// m -A a +iO \ 
^ sinh JV (A a -iC/2) r ii 1 sinh( Mm -A a -iC) y l ' 

A M ({n}, {A}) = /( 7' C/2 ^ fl sinh (^ - </ 2 )' (E.35) 
sinni — iC) ± \ 

v ,y m=l 

, n M • w\ \ M-l 

n /r 1 rxi\ 1 llm-1 Sinh(A 6 - Hm) -pr . 

sinn(A 6 + iQ/2) n m =i,^ 6 smh (Ab - A m ) m=1 ^ b 



,(E.36) 



sinh(A b - i(/2)f(\ b - iQ sinh(A b + 3</2)/(A b + iQ 
sinh(A Q - A b + iQ sinh(A a - A b - iQ 

B Mb (M, {A}) = (/(A b + iQ f(X b iQ) JlF= lSinh ^'^\ 

M-l 

x Yl sinh(A m -A b -<), (E.37) 

m=l,^6 
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- - J^-lft) n iigi n + -C/2XE.38, 
%m(Mi!M) - n S^rfi n + < E - M > 

By setting /(z) = sinh(z — i(/2) we obtain 

M 

det(S ab A a + B ab ) = Y[Ai det (<U + G J , (E.40) 



i=l 



where the (M — 1) x (M — 1) matrix G ab is given by the following elements 

q- = 1 IlmLi sinh(A b - fi m ) IlmJi sinh(A m - \ b - i() 

ab I sinh (A 6 + </2)l 2 n™=i,^6 sinh ( A " - A ») n„=i sinh( Mm - A 6 - iC) 

x L sinh jv (A fc + </2) yr sinh(/x m - A b + iQ \ 
X ^ sinh^Ab - iC/2) sinh(/i m - A* - <) J 

x / sinh(A b - iC/2) sinh(A b - 3»C/2) g \ 
V sinh(A a -X b -iQ ' )' 

The final expression reads 

j M M-l 

d M ^ = smh^^vzn n n sinh (^« - ^ - o ^ + g ^) 

^y-r 1 / sinh w (A t + </2) -pr sinh(^ m - A^ + iQ \ 
11 1 sinh N ( A, - i(/2) sinh( Mm - A, — *C) / 

x nfiT 1 I sinh(A t + »C/2)| 2 n,>, sinMMj - Mi ) Yjf^ 1 sinh(A J - - AQ 
" I ^Mrt + </2)| 2 ~ IL=i 11*17 1 sinh( Mi - A,) 

Appendix F. Thermodynamics limits for the xxz spin chain 



x 

i=l 



(E.42) 



In this Appendix we consider the thermodynamic limit of various expressions that are 
part of the form factors of the operators S* and S~ of the xxz spin chain. 

Appendix F.l. Calculation of Mi 

Let us repeat one more time the formula for Mi 

Mi= f[ sinh( Mt - A 3 - iQ ^ (F.l) 

i sinh 1/2 (^ - - iQ sinh 1/2 (Ai - Xj - iQ 

The crucial point in the calculation of double products like Mi and Mi is to divide 
them into 3 different parts as follows 

_ _ TT" sinh( M , - Xj - iQ 

i j=l sinh ' [fii - Hj - iQ) smh ' (A; - Xj - iQ 

p2 = TT TT » sinh (/V - Mm - <) Sinh(^+ - X m - <) 

2 f=i i=i sinh (M 4 + - Mm - <) sinh(/ir - A m - iQ 
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™ sinh(^+ - fjq -iQ 

P3 = 11 K) s^ /2 (p7 7j <0 ' (F ' 4) 

In terms of these Mi is simply 

M l =p lP2 p3. (F.5) 

In p\ we can take the thermodynamic limit of both products, in p 2 only of the inner 
product, and ps is a discrete part which depends on the details of the excited state. 
We start with pi , which gives 



M 



sinh(/Xj — \j — iQ 



i sinh 1/2 (^ - - iQ sinh 1/2 (Ai - Xj - iQ 



(F.6) 



sinh ( \i — Xj — iQ — 



FnVh) 
Np N (\i) 



M 

Ui -nh 1 / 2 (A, - Xj - <C - + sinh 1 / 2 ( Aj - A, - <) 

™ cosh (wivh) 1 - coth(A t - A, - iQ tanh (jjgfo) 

11 , 1/2 / F„fA,1 fVfA,) \ / „»v\\l/2' 



'5=1 ^h 1 / 2 (JM^ - ( x _ coth(At _ A . _ <} tanh _ 

In the following calculations we will again skip some expression that contribute only 
to the phase of the form factor. The first product gives 

F N (\i) 



cosh 



Np N (Xi) 



cxpQy dXd^F(X)F^yj +0(1/N). (F.7) 



M 

.11 ~ ,1/2 / F N (M) F N (X 3 ) \ 

JJ= icosn y NpN{Xi ) Np N (X j )J 
The second product gives 

l-coth(A i -A j -<)tanh( 1 ^ y ) 



M 

n 



exp 



x exp 



x exp 



= exp 



(l - coth(A 4 - Xj iQ tanh - gjM)) 

1 Fn(Xi) 



1/2 



(F.8) 



2 ^ Np N (Xi) 



(coth(A^ — Xj — i() — c.c) 



1,3 



FN(Xi)F N (Xj) 



coth 2 (A, - Aj -<) +C(l/iV) 



2 N 2 p N (Xi)pN{Xj) 

dndXF N (X)F N {n) (coth 2 (A -ii- iQ) + i<f> ) + 0(1/N), 



1 /■« 
2 



where in we gathered all the purely phase factors. They do not contribute to the 
final answer and can be dropped. Note that despite factor i the integral is a real 
number. Altogether we obtain for p\ the following answer 

\ Pl \ = exp (-\ f +0(1/AT). (F.9) 



tanh (A — /i — i£) 
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We move on now to p 2 



M 



P2 =nn - 



sinh(/x t - n m - i() smh(fj,f - X m - i() 
1 sinh(/x+ - p m - K) sinh(/ir - A m - <) 



82 



(F.10) 



™ M sinh (it- - A m - < 

mi 



Fw(A m ) 
Np N (X m ) 



sinh(yu+ - A m - i() 



i—1 m— 1 



n M 

nn 



sinh(^ - A m - <) sinh (^if - 
1 + coth(^ - A m - <) sinh ( 



, Fjy(A m ) 



Fiv(A m ) 
WpjvfA™) 



+ 0(1/N) 



i\ m= i 1 + coth( Ai + - A m - iC) sinh ( 1 g^y) 
™ / /•<? 

Ill cxp / dAF ( A ) (coth(/i,r — A — »C) — coth(^+ — A — i()) ) . 

i=l ^■'-9 



Which gives 



P2 



V J -9 tanh(,u--A-iC) / 
1 \ J -q tanh(/i+-A-»C) / 



« exp 

n 



The final answer for Mi is then given by the following expression 

2 



\Mi\ = cxp 



9 d,dx F ^ F ^ 

q <.,,.!.-< 



tanh (A — fi — iQ 



(F.ll) 



(F.12) 



„ expff dX h /j A \ 



,i exp 



(A 



dA 



><n 



£w 

tanh(/j.+ — A— i£) / 

sinh(/x+ - M7 - O 



sinh 1 / 2 ^ - M + - »C) sinh 1 / 2 ^ - M T - <) 
A similar formula holds for Mi with f+(A) replacing F(A). 



■0(1/N). 



Appendix F. 2. Calculation of M 2 

Let us recall once again the expression we are going to consider in this section 

sinh 1 / 2 ^ — Hi) sinh 1//2 (A., — Aj) 



M 2 =n ; 



sinh(yUi — Aj) 

The double product can be written as a product of 3 different expressions 

n M „sinh 1/2 (^ - Hi) sinh 1/2 (A i - Aj) 




sinh(/ij — Aj) 



n M 



*=nn 



sinh(A t - fi h ) sinh(/Xi - 



t=i Li sinh (Mi - V h ) sinh(A, - fj+) ' 
sinh 1 / 2 ^ - M +) sinh V2 (/U 7- M 7) 



sinh(^+ - /x • ) 



(F.13) 

(F.14) 
(F.15) 
(F.16) 
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which we obtained from rewriting the expression under the double product as 

„sinh 2 (/i i - nf) 



83 



Y[ sinh(^- - m) = Y[" sinh(>j - Mi) IJ Yi 
J^j sinh(/i^ — fij) sinh(/x~ — fj,J) 



A I 



sinh 2 (^+ - fij 

M 



n M 



! [ sinh( M i — Aj) = []" 8inh (Mi - A,) If [] 



Let's focus on the first expression. We start with 
l,/ sinh 1 / 2 ^^ — (tx^) sinh 1 ^ 2 (A J — A.,) 



»=i j'=i#« 



sinh 2 (^ ) 



, sinh(^+ - Aj) 
sinh(^r - Aj) 



P1 =n 



sinh(u ■ - Aj) 



" ^V 2 - ggj (l - cotMA, - A,) tanh (j| 



)_ _ F(A,) 
JVp(Ai) 



1/2 



1 — CO 



th(A i -A,)tanh( 5 g^ y )) 



exp 



d\F(X) ] 



n 



(l-coth(A J -A l )tanh( 1 g^ y 



1/2 



1 — CO 



th(A i -A i )tanh( 1 g^)) 



where we used that 



M cosh V2 

cosh( — — 



)_ _ F(\j) 
) Np(\ t ) 



exp 



cZAF(A) J 



0{l/N). 



(F.17) 



(F.18) 



Calculation of the second product is more difhcult. Let us use the following notation 

1/2 

x \ i ( .,,.1, I ' 1 1 i 



m ( 1 - cotMA, - Ai) tanh ( ^ - ^ 



(l-coth(A,-A l )tanh( 1 g^ T )) 



(F.19) 



and consider different regions. We follow the same logic as in the case of the Bose 
gas. The names of the regions are also the same. The details of the computation are 
outlined in section FOl 



Appendix F.2.1. Region I In the first region, where we can expand the product 
logarithmically, we obtain 

l°S*> = i £ coth^-A,)-^"^ 



2 y 
1 ri 



q+v*{-q) 



dX 



dXF 2 (X) 



n 2 p(xm^ 

rq-v*(q 
d/1 + I dX 



d\x coth 2 (Aj - X, t )F(X)F(p). 



A+i/*(A) 



A-i/*(A) 



d\i coth 2 (A — (it) 
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r q-v*(q) r q 



dXF 2 (X) f dn coth 2 (A - 

i rq r\-v*{\) 

- dX dfi coth 2 (A - ft)(F{\)- F(fi)) 2 

J — Q+V ( — a) J — q 



-q+v*(-q) J-q 

i rq-v*(q) rq 

-- dX rf/xcoth (A — /i) (-F(A) - F(fj,)) . (F.20) 

4 J-g J\+v*(\) 

The last two integrals are actually cut-off independent, one can extend limits of 
integration without introducing any error nor divergence, therefore 

i rq r \-v*(\) 



2 



r dXFHx) f 

J -q+v*{-q) J-q 



logT/ = -/ dXF\X) djj, coth. z {X- ijl) 



1 rq-v*(q) rq 
+ - / dXF 2 {X) / d M coth 2 (A- J u) 

2 J-q J\+v*(\) 

-j f 9 dX f 9 ^coth 2 (A - fj) (F(X) - F(n)) 2 . (F.21) 
J —q J—q 

The first two integrals are cut-off dependent and should simplify once we combine 
answers from different regions. Now let us move to region II. 

Appendix F.2.2. Region II In this region we can make the following approximation 
where 



x > -»' = is 6 - *> - svfe {i - kf + OWN ' % (F - 22) 



C > = -|M. (F.23) 



First one observes that 

F(Xj) F(Xi) _ d fF(Xj) 



|A=A,(A J -A i ) = C;(A j -A i ),(F.24) 



p{X 3 ) p(Xi) dX V p{^j) 

where we have defined Cj. We only need to consider terms up to 0(1/N) because 
they are the only ones contributing in the thermodynamic limit. Then we have 

Tii= j-j tanh 1 / 2 (A, - AQ (tanh(A J - -AQ- tanh (g^ - g^)) ^ 
i^-e/j (tanh(Aj - A*) - tanh (w^)) ) 

tanh(A, - Aj) tanh ( A J ~ A ~ tanh (jgAj) ~ Wx6)) 

K i^ jeII tanh(Aj - A 4 ) - tanh (ig^jy) tanh(Aj - A,) + tanh (^jj 



For the first product we obtain 

tanh(Aj — Aj) 
tanh(A, - A,) - tanh ~ j-i- 

1 Ki 



JC 6 



0(l/iV 2 ). 
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For the second 

tanh(A, - A*) - tanh ~ wMl) _ j-i 



tanh(Aj - A;) + tanh (jf^y) 



. _ . F(Aj ) 



+ 0(1/N 2 ) (F.27) 



l- 



_J 2_ 



Therefore 



T »= II ZWTT~ (F.28) 



We obtain the same expression as for the Bose gas (see Eq. (II 19|) , therefore as a final 
answer in this region we have (Eq. (|126[1 ) 



x exp -— / dAp(A)F 2 (A) - — / dAp(A)F 2 (A) 



r(i - F( 9 ))r(i + F( g ))r(i + F(-«))r(i - f(-<?)) 



i*/2 

. (F.29) 



Appendix F.2.3. Region III In region III we again obtain the same answer as for the 
Bose gas (Eq. ([125]) . Eq. (fT2Ul) ). hence 

T 7W = [G(l + F(g))G(l - F(q))G(l + F(-q))G(l - F{-q)f' 2 

x (n* y F2{q) ^ q) [r(l - F(g))r(l + F(«))r(l + F(-«))r(l - F{-q))f' 2 . (F.30) 

Appendix F.2-4- Cutoff dependent part In the end, let us collect all the cut-off 
dependent parts of Tj, T/j and T///, altogether they give 

log(T cutoff ) = -\{F 2 {q) + F 2 (-g)) log(n*) 

+ - / d\F 2 {\)[ ^ co th 2 (A- M )-— - 

1 rH-v*{q) / rq i \ 

+- / dXF 2 (X) I / d M coth 2 (A- M )-— - . (F.31) 

Let us start with one of the inner integrals 

/A-i/*(A) i o 
d/jcoth 2 (A - /i) -— = A - -i/* (A) +q- coth(A + g) + O ({v*(X)f) 
-q V (A) 

(F.32) 

and then 

dA F 2 (A) / d/j coth 2 ( A - /i) — 

-<;+*>"•(-<?) \J-q 



q+v*(-q) 

dAF 2 (A)coth(A + g) + / dXF 2 (X)(q + A) + 0{v* (A)) 

-?+"*(-«) 
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+ f 9 dXF 2 (X)(q + X) + <D{v*(X)). (F.33) 

J-q 

Similar calculation for the second integral gives 

/q-v*(q) ( r-q i \ 

d\F 2 (\) I / ^coth 2 (A-/x)-- 7 TT 

/q-v*(q) r-q 
dAF 2 (A)coth(A- g ) + / F 2 (X)(q- \) + 0(v*(X)) 
-q J —q 

+ f 9 d\F 2 (\)(q - A) + 0{v*(X)). (F.34) 

J-q 

Finally for the cut-off dependent part we get (p(q) = p{—q)) 

log(T cutoff ) = - l -{F 2 {q) +F 2 (-q))\og(qNp(q)) +q J" d\F 2 (X) 



+I Pj _ C dx 1 F2 ^ Ip C dx ( F 35) 

2 + J_ 1 tanh(g(z - 1)) 2 J_ x tanh(g(a; + 1)) ' v ' ; 

Gathering together results from all the regions we obtain 

(F 2 { q ) + F 2 (- q )) 



M 

X 
3 



Pl = (qNp(q)) — — (G(l + F(q))G(l - F(q))G(l + F(-q))G(l - F(-q))) 1/2 

x exp Q j" d\(l - ,F(A)mt(rf(A))) (V(A) - F( ^' ) W )) 

'p V fc «fV> _ 'p. /' ix , (F.36) 

2 tanh(g(a; - 1)) 2 y_! tauh(q(x + 1)) J ' v ; 



-9 

x exp ( ^-P. 



where we used that 



q J" d\F 2 (X) = 1 ^ d\dp (F(X) - F{p)f + l - ^ dXF(X^) (F.37) 

Appendix F.2.5. Calculation of p2 The calculation of pi is divided into 2 steps. First 
we will calculate the hole terms, thereafter the particles terms. We start with 

= IJ I| l-coth(A i - M ,)-^- x(l + CJ(l/JV)). (F.38) 

There are again two qualitatively different situations, first when holes are near the 
Fermi edges, second when they are located deep inside the Fermi sea. Without making 
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any mistake we can consider a single hole and then simply generalise to many holes by 
multiplying single contributions. Divide now the product into two parts, first when 
the hole rapidity is well separated from other rapidities, and second when the hole 
rapidity is close to the other rapidity. We obtain 

TLe = expf £ COth %- X 'f iX A + 0(l/N), (F.39) 

\\j-h\>n* P( z) J 

TlL= II l ~\ + _ F h iXh) +0(l/N), (F.40) 

\j-h\<n* 

T hole =Tl ole xTH le . (F.41) 

Case I : Hole is near the quasi-Fermi point 

We obtain a similar expression as for the Bose gas (Eq. (|138[) ) 

T h,edge _ , „,«tf (±«) sin (^F(g)) T (TQh ± F(±q)) 

hole -wm) nF(q) r(Tg ± } 

xexpfp ± f dx g^L V (F.42) 
V 7-1 tanh(g(a;Tl))/ 

Number is a quantum number that parametrizes the position of the hole and is 
defined in the same way as for the Bose gas. 

Case II : Hole is deep inside the distribution 

The same situation here (Eq. ([13610 

Tttr = ^^^^^ exp [ P f dx f( qX \._^ | • (F.43) 



* F (Ph) I J-i tanh (q (x - ^ 



Particle term 

The computation of the particle part goes along the same lines and gives (compare 
with Eq. ([TIP]) . Eq. (fTTO l 



"^- p ''> J(^l)) )- (F - 44) 

1 tanh I q ( x - — 



rph,far 

1 particle ~ eX P 

Here again p^ parametrizes the position of a particle in the same way as for the Bose 
gas. 

Appendix F. 2. 6. Final answer Let us summarise now the calculation of M2 . In order 
to follow the same notation as in the case of the Bose gas we redefine term p^ as T cross 
(see Eq. fug) ) 

nsinh 1 / 2 ^ - nf) sinh 1 / 2 (/x~ - 
3 . /? + — • (F.46) 
m smh(/i+ - (i. ) 
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The final answer for M 2 is 

M 2 = (qNp{q))- {F2{q)+ ^ ( ~ q)) (G(l + F{q))G{\ - F(q))G(l + F(-q))G(l - F(-q))) 1/2 

\ J i r jw ../f(a)-f(/x) 



x 



n(^H4/>( 

cxp Q y 9 dA (1 - ttF(A) cot(7rF(A))) (V(A) - 



x 



sinh(A — fj,) 

P(A) 



! tanh(g(x - 1)) 2 ,/_! tanh(q(x + 1)) _ 

n 

X IT (^/iote X T< plrUcle ) x T cross . (F.47) 
i=i 

For the form factor of S~ operator calculation of M 2 follows the same lines and 
the final answer is the same up to a change in shift function, F + (X) instead of F(X). 

Appendix F.3. Fredholm determinants 

The last crucial step that wc need to perform is take the thermodynamic limit of 
matrix elements G z ab and write its determinant in form of a Fredholm determinant. 
The expression for G z ab reads 

sinh«) IIm=i sinh(A 6 - p m ) yr sinh(A m - A b - <) 
ab sinh(A fe + iC/2) n ^ =1 ^ s inh(A 6 -A ro )ii 1 sinh( Mro -A 6 -iC) 1 ' ) 

x / | sinh jv (A a +2C/2)n^iSinh( Mm -A b + zC) y 1 
^ sinh^Ai, - tC/2) n„=i sinh( Mm - A 6 - »C) / 

/ sinh(A b + 3</2) sinh(A 6 - iQ/2) ji sinh(A m - </2) 2sinh«) \ 
ysinh(A a - Ab - <) sinh(A a - A 6 + <) sinh(/i ro - </2) sinh(A a - i(/2) J 

We start by considering the first product 

IIm=l Smll ( A a - Mm) _ , „ A TT sillh ( A a ~ Mm) 

. 8inh(A „ _ „„, ft n ^_^a) 

Ai smh(A a - /i, ) A ^ smh(A a - A m 

In order to calculate the thermodynamic limit of the last expression we need to be 
careful and split calculations into two regions. Let's define the cutoff distance n* and 
associate to it two rapidities Al and A_r corresponding to rapidities with indices a — n* 
and a + n* respectively. Then we have 

" smh(A a - A m ) ax \^ Np(\ m )J 

= n (i + ^ - a„,,m«) n (, + „ ti ,A„ - a„,,^i) 

|m— a|>n* |m— a|<n* 
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= exp ( / d\F(X) coth(A a - A) + / dXF(X) coth(A a - A) 



TT ( a - m + FNiXm) )+o(i/N) 



x 

\m— a\ <n* 



r(l + n* + F N (X a ))T(l + n* - F N {X a )) 



= exp (p£ dXF(X) coth(A Q - A) ; p2(i + ^ + _ 
= exp (p J" dXF(X) coth(A Q - A)) + 0{\/N). (F.50) 



Therefore 
Ilm=i sinh(A a - fi m ) _ sin(7rF(A a )) 



exp ^P J dXF(X) coth(A a - A) 



n^=i # „smh(A a -A m ) ^Vp(Aa) 

x ^ sinh(A a -^) 
l= i smh(A a - /ij ) 

Now we move on to the next product. Here the thermodynamic limit is straightforward 
and gives 

yr sinh(A m - A a - <) _ -A- sinh(^~ - A a - i() 
11 sinh(/i m - A a - iC) Al sm h(>+ - A„ - <) 

x exp dXF(X) coth(A - A a - <)^j + <D(1/N). (F.52) 
The next expression can be written using Mo defined earlier Eq. (|225[) . and reads 

(i + Moc^))- 1 = (i + sinh :^ + ^n^ 1 ^^A a+ < ) \ - 1 

y smh (A a - i(/2) n m =i sinh(/i m - A a - if) / 

Finally, the last product gives 

sinh(A m - i(/2) 



^ S inh(/i m + iC/2) 

A sinh(u7 + iC/2) /f 9 F(X)dX \ , , , , 

= TT r-T ^rrexp / ; — r +0 (IN . F.54 

sinh( M + + if/2) P U-,tanh(A + </2)y w ' V ; 

Therefore, the matrix elements G ab take form 

GZ _ 1 sinh«) -A- sinh(/^~ - A a - Q sinh(A a - 

ab ~ 2irNp(X a ) sinh(A b + </2) 11 sinh(^+ - A a - <) sinh(A a - n~) 

x exp ( P / 9 dX -ffi^ + / 9 dA f ^ — + i7rP(A Q ; 



tanh(A a — A) J_ q tanh(A — A a — i() 



sinner - i(/2) ff q F{X)dX \ -2sinh«) 



JJ sinh(^ - iC/2) 6XP V./ « tanh(A-tc - ' siiih(A„ - k,/2) 

ih(A 6 + 3</2) sinh(A fe - iC/2) ' 
sinh(A Q - A 6 - iQ sinh(A a - A h + i() / 



sinh(A + 3</2) sinh(Aft - </2) 
+ :ZX7^ x 77T - ~ZTT7"\ x I 7F\ • ( F - 55 ) 
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Now we can define an operator G that acts on (— <?, q) and is given by 

fr sinh(/i- - i'C/2) / f g F(X)d\ \ -2sinh«) 



A,, x , x [ A smhQy-jC/2) //« F(X)d\ \ -2 S inh(iQ 
(M ' "> = ^ (J ! sinh(Mi+ _ iC/2) ex P (X, tanh(A-</2)J sinh^ - </2) ' 



sinh(i/ + 3iC/2) sinh(f - i(/2) 



sinh(/i — v — iQ sinh(/i ~ v + iQ 
where 

1 sinh(i£) -p-r sinh(//,~ — f — iQ sinh(V — /i^ 1 



(F.56) 



_ i aumyt^ T-r 



27T sinh(V + «C/2) sinh(^+ — v — iQ sinh(V — ) 

x exp (p [ 9 dX F / A) - + I" dX — ^ — + inF(u)] . (F.57) 

V tanh(f - A) J_ q tanh(A - i/ - i£) K ' ) v y 

Finally 

det(5 o6 + G z ab ) = dot (l + G z ) + 0(1/JV), (F.58) 
where it is understood that the determinant on the r.h.s is a Fredholm determinant. 

Appendix F.3.1. Determinant in S~ form factor The matrix appearing in the 
determinant in the S~ form factor is given by Eq. (|212j) 

q- = 1 Um=i smh(X b - fj, m ) Um=i sinh(A m - X b - iQ 

ab I ^h(X b + </2)P n ^ =i ^ sinh(Ah - X m ) sinh(^ m - A 6 - iQ 

x f 1+ sinh jv (A fc + </2) -pj sinh(/i m - A a + <) \ 
X ^ smh N (X b ~ i(/2) sinh(/x m A a i0 J 
sinh(A b - </2) sinh(A 6 - 3</2) 



c.c. . (F.59) 
smh(A a - X b - iQ J 

The thermodynamic limit here is similar to the case studied above. The only difference 
is the extra rapidity (UjLj appearing now and then. In the thermodynamic limit we 
have 

det (S ab + G~ b ) = dot (1 + H-(jjt, v)) + 0(1/N). (F.60) 
The operator H~ (ji, v) acts again on {—q, q) and is given by 

H = b z/) — — c.c , (F.61) 

\ smn(/i — v — iQ) J 

where 

b -, v s _ 1 A sinh(^~ -v-iQ yf sinh(^+ - v) 

{) 2tt| sinh(zy + </2)| 2 11 sinh(^ - v) smh{^f - v - iQ 

x exp ( [" dX ^ +(A) - - P [" dX f+ (A) , + iirF + {X)) . (F.62) 

\J- q tanh(A - v - iQ J- q tanh(A - v) J 
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Figure 6. Log-log plot of A+(k) (black, dashed), A~(k) (orange), A + (k) (blue, 
dashed) and A~ (k) (red) as functions of k/kp for K = 1.81. 
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Figure 7. (Left panel) Plots of Di and D A (c.f Eq. I|277|) ) as a function of A for 
two values of s z (average magnetisation per site s z = 0.5 — M/N). The prefactor 
D\ approaches the free fermions limit as A — i- 0. (Right panel) Plots of Di as 
function of s z for a few values of A. The results agree with the numerical results 
of [49]. 
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Figure 8. (Left panel) Plot of Eo and E\ (c.f Eq. (|278jl) as a function of A for 
different values of s z . (Right panel) Plot of Eo as a function of s z for different 
values of A. Numerical data (crosses) from Ref. |49| . 



